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1.0  INTRODUCTION 


The  simulation  facility  currently  operational  at  the  Guidance  and 
Control  Directorate,  Missile  Research  Development  and  Engineering 
Laboratory,  U.  S.  Army  Missile  Command  consists  of  a hybrid  digital/ 
analog  computer  system  Interconnected  to  a pair  of  target  simulators  and 
to  a three  degree  of  freedom  angular  motion  simulator.  The  hybrid 
computer  comprises  a medium  scale  digital  computer,  as  the  controlling 
eluent  of  the  facility.  Interfaced  to  a pair  of  analog  computers;  target 
simulators  permit  the  simulation  of  laser  designated  targets  and  optically 
tracked  targets. 

By  placing  guidance  and  control  hardware  components  on  the  angular 
notion  simulator  which  Is  driven  to  simulate  the  closed  loop  missile 
angular  response,  and  arranging  for  the  target  simulator  to  generate  the 
appropriate  cloSed  loop  line  of  sight  angles  and  angular  rates,  the  per- 
formree  of  the  actual  flight  hardware  can  be  studied  and  evaluated.  This 
requires,  of  course,  a real-time  solution  of  the  missile  flight  equations 
to  be  provided  by  the  hybrid  computer. 

Descriptions  of  the  simulation  facility  hardware  components,  their 
performance  characteristics  and  methods  of  utilization  are  scattered  In  a 
diversity  of  reference  and  Instruction  manuals  which  do  not  lend  themselves 
to  yielding  a rapid  overall  grasp  of  the  capability  and  operation  of  the 
facll'ty.  The  Intent  of  this  report  Is  to  remedy  this  situation  by  pro- 
viding a summary  of  the  list  of  components  which  comprise  the  “facl l ity,  a 
brief  description  of  tlielr  characteristics  and  methods  of  use  and  any 
other  salient  points  of  particular  interest  to  the  potential  user,  it 
should  be  noted  that  the  hybrid  computer  software  components  are  considered 
essential  parts  of  the  simulation  facility. 

A large  part  of  the  Information  contained  In  Sections  2 and  3 Is 
available  elsewhere,  and  references  are  given  wherever  appropriate. 

However,  some  of  the  information  Is  not  generally  available  In  written 
form. 


Finally,  an  example  hybrid  simulation  program  Is  given  which  con- 
sists of  a missile  Intercepting  a target  with  both  restricted  to  planar 
motion  (3  degrees  of  freedom).  A discussion  of  the  method  of  solution 
of  the  problem  Is  given  to  Illustrate  the  use  of  the  hybrid  computer  and 
to  Indicate  a general  approach  to  setting  up  a real-time  simulation. 


2.0  COMPUTER  EQUIPMENT 

2.1  Digital 


The  digital  section  of  the  hybrid  computer  facility  consists  of  a 
Xerox  Sigma  5 digital  computer  system.  The  system  configuration  Is  shown 
In  Figure  2.1,  and  a list  of  the  system  components,  Including  Xerox  model 
numbers  and  brief  performance  specifications,  Is  given  In  Table  I.  The 
table  does  not  Include  details  of  the  analog  computer  Interfaces;  these 
are  given  In  Section  2.3  below. 

2.1.1  Central  Processor  Unit 

The  operational  environment  within  the  CPU  Is  controlled  by  a 64 
bit  Program  Status  Word  (PSW),  which  allows  rapid  switching  of  the  environ- 
ment In  response  to  external  stimuli  or  under  Internal  program  control. 

The  central  processor  Instruction  set  Includes  control  logic  for  execution 
of  word  (32  bits),  double-word  (64  bits),  half-word  (16  bits),  and  byte 
(8  bits)  oriented  processor  commands.  Additionally,  a f‘.  atlng  point  unit 
handles  arithmetic  on  single  precision  (32  bits)  and  double-precision  (64 
bits)  operands.  Floating  point  format  comprises  a biased  8 bit  exponent, 
with  base  l 6,  at  the  most  significant  end  of  the  word  or  double-word,  and 
a 24  bit  or  56  bit  mantissa. 

General  purpose  registers,  which  are  used  for  all  arithmetic,  logic 
and  Indexing  operations,  consist  of  blocks  of  16  32-blt  words/  The  current 
configuration  contains  4 register  blocks  (numbered  0,  1,2,  3)  with  a 
capability  of  expanding  to  16.  Register  block  usage  Is  controlled  by  the 
Register  Pointer  bits  In  the  PSW. 

All  computer  Instructions  occupy  32  bits  v.Jth  a 17  bit  address  field 
which  allows  direct  access  to  the  maximum  memory  address  range  of  0 to 
131*071  words.  Addressing  modes  Include  direct,  Indexed  (using  1 Index 
register).  Indirect  (1  level).  Immediate,  and  combined  Indexed  and 
Indirect  (post-indexing). 

The  system  contains  two  real-time  clocks  with  basic  clock  rates  of 
8000  per  second  and  500  r second  respectively  (l.e,  resolutions  of  .125 
and  2 milliseconds),  ftor1  Interrupts  may  be  generated  such  that  a user- 
chosen  multiple  of  the  r'slc  rates  may  be  maintained.  Details  of  the 
methods  of  use  and  programming  at  the  CPU-hardware  level  are  given  In 
Reference  I . 


2*1.2  Input/Output 

As  Indicated  In  Flqure  " 1 Input/output  operates  through  three  paths 
connecting  magnetic  core  memor*  with  the  1/0  devices. 


FIGURE  2.1  Sigma  5 System  Configuration 
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TABLE  I 

Elements  of  the  Sigma  5 Digital  Computer  System 


Description 

Xerox 

Model 

Number 

Performance 

Characteristics 

, " — — - » 

Ml OP  Channel  ! 
Assignment 

Sigma  5 CPU 

8202 

See  Ref.  1 

- 

NIOP 

8273 

400  KB/sec. 

i 

i 

Core  Memory 

Various 

.85  usee  cycle 

i i 

i 

Card  Reader  1 

7122 

400  cpra 

3 ! 

1 

Card  Reader  2 

752i 

200  cpm 

4 1 

Line  Printer  l 

1 

7*41 

1100  1pm 

2 

Line  Printer  2 

7440 

600  1pm 

1 

1 

5 ! 

i 

; Keyboard/Printer 
, (Teletype  ASR-35) 

7020 

10  char/sec 

i 

1 

i 

RAD  Controller 

7201 

188  KB/sec 

- Fx  ! 

RAD  File  0 

7204 

3 Megabytes 

F0  j 

RAD  File  1 

7204 

3 Megabytes 

FI  1 

i 

Magnetic  Tape 
Controller 

7320 

60  KB/sec 

Bx  ! 

1 

Magnetic  Tape  Mve 

7322 

75  lps 

80  j 

Graph  Plotter 
Controller 

7534 

.01  Inch  >nc 

7 j 

1 

Communications 
Control ler 

7611 

1800  Baud 

F 

CRT  Display 

4002A* 

600,  1800,  2400 
or  §600  Baud  1 

- 

•Tektronix  Inc.  Hods!  Number 


TABLE  II 


DIO  Address  Characters  (8  Bits) 


1 

1 ■ 

Address  Character  . 

1 Device  Description 

(Upper  8 bits)  | 

Cl-510  (Interface  to 
Cf-5000  Analog) 

CI-500  Interface  and 
DttA  Control  Signals 

Cr'innur."! cations  Controller* 
( . m litter  Output) 


* The  controller  Is  numbered  1 so  that  all  address  fields 
are  of  the  form  1 30 lx*  where  x specifies  the  required 
function. 
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Data  transfer  operations  via  the  DIO  (Direct  Inpu'./Output)  takes 
place  through  CPU  registers  and  occupy  the  CPU  for  the  duration  of  the 
transfer.  Each  word  transferred  requires  a 16  bit  address  field  In  the 
computer  Instruction  to  identify  the  I/O  operation  and  the  device  involved. 
Table  II  shows  the  most  significant  byte  of  the  address  field  for  the  hard- 
ware currently  attached. 

The  H I OP  (Multiplexer  Input/Output  Processor)  transfers  data  In 
parallel  with  CPU  operation  and  occupies  the  CPU  only  during  the  time 
required  to  transfer  control  signals  to  or  from  the  CPU.  The  MIOP  can 
handle  up  to  32  I/O  channels  each  operating  simultaneously,  provided  the 
overall  data  transfer  rate  does  not  exceed  the  MIOP  maximum  rate  band- 
width of  about  b 00  kl lobytes /second . If  the  maximum  rate  is  exceeded  at 
any  time  then  a 'data  overrun'  error  condition  Is  signalled  for  some  or 
all  of  the  peripheral  devices  involved. 

The  MIOP  currently  contains  16  data  channels,  of  which  9 are  assigned 
to  peripheral  device  controllers.  Two  further  groups  of  8 channels  (Xerox 
model  number  8276)  may  be  added.  The  Graph  Plotter  controller  shown  In 
Figure  2.1  Is  Intended  to  drive  a 12  Inch  Calcomp  plotter;  however,  the 
plotter  Is  not  currently  attached  to  the  system. 

Data  transfer  to  and  from  the  communications  controller  (Xerox 
model  7611)  uses  both  the  DIO  and  MIOP  data  paths.  Input  is  via  the  MIOP 
and  output  Is  through  the  DIO  on  a character-by-character  basis  with  an 
external  Interrupt  generated  following  the  transmission  or  receipt  of 
each  character.  Note  that  the  communications  controller  does  not  use  the. 
Internal  I/O  Interrupt  as  do  all  other  devices  which  use  the  MIOP.  Data 
transmission  speed  for  the  single  communication  channel  currently  used  Is 
set  at  9600  Baud.  The  7611  has  the  capability  of  expansion  to  6b  channels. 

The  third  I/O  path  consists  of  a Direct  Memory  Access  channel  which 
Is  used  exclusively  to  connect  the  ADC's  and  DAC's  In  the  Interface  to 
the  Cl -500  Analog  computer  to  core  memory.  This  channel  operates  In 
parallel  (asynchronously)  with  CPU  operation  but  Is  controlled  by  signals 
transmitted  via  the  DIO. 


The 

conversion 
of  any  I/O 
Interrupt, 
accessible 
and  DID  to 


DMA  channel  operates  at  high  speed  (limited  by  ADC  and  OAC 
rates),  transferring  data  In  32  bit  parallel  format.  Completion 
operation  may  be  signalled  to  the  CPU  by  use  of  an  externa! 

It  should  be  noted  that  the  CI-500  ADC's  and  DAC's  are  also 
via  ths  DIO  channel.  However,  It  Is  not  possible  for  the  DMA 
access  DAC's  or  ADC's  simultaneously. 
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2.1.3  Magnetic  Core  Storage 

Magnetic  core  storage  consists  of  49,152  (48K  where  K * 1024)  words 
of  memory,  divided  Into  3 banks  of  16,384  wordr.  each.  Banks  are  grouped 
Into  2 memory  units,  each  of  which  have  separate  access  ports  for  the  CPU, 
MI0P  and  DMA.  Address  Interlacing  between  banks  and  units  Is  used  to 
achieve  overlap  In  accessing  particular  memory  locations. 

Memory  locations  0 through  15  are  not  available  since  any  reference 
to  these  addresses  In  a computer  Instruction  refers  to  the  general  register 
of  that  number  In  the  block  selected  by  the  Register  Pointer. 

2.2  Analog 

The  analog  section  of  the  hybrid  computer  facility  Includes  two 
cnclog  computer^  - Astrodata  Comcor  C l -5000  and  Astrodata  Comcor  CI-500 
machines.  These  machines  are  solid-state,  general-purpose  analog  computers 
with  operating  ranges  of  +100  volts.  They  each  Include  a good  range  of 
linear  and  non-linear  anaTog  elements  plus  powerful  logic  control  sections, 
A brief  dtscussion  of  the  analog  and  logic  components  In  both  machines  Is 
given  below. 


2.2.1  Analog  Components 

Operational  Amplifier 

The  operational  amplifier  may  be  used  as  a combination  summer/ 
Integrator  with  Input  gains  of  10  or  1 available.  Additional  gain  as  an 
Integrator  may  be  obtained  by  choice  of  capacitor.  Gain  due  to  a capacitor 
Is  equal  to  vC  (since  Input  resistors  are  1 megohm)  where  C Is  given  In 
micro-farads.  Capacitors  available  are  10  uf,  1 uf,  .1  uf,  .01  uf  and 
.001  uf.  Only  odd  numbered  amplifiers  may  be  used  as  Integrators. 

Capacitor  selection  Is  determined  from  the  logic  board  as  shown  In 
Figure  2.2. 

Referring  to  Figure  2.2,  amplifier  003  Is  connected  as  an  Integrator 
using  electronic  switching  (S^  to  A003)  In  the  real-time  mode  (RT  to  .1) 
with  gain  of  10  (*/C).  Reset  and  hold  conditions  are  determined  from  the 
keyboard  (RR  to  R,  HH  to  H).  In  order  to  use  relay  switching  remove  (S^ 
to  A003).  For  the  timescale  mode  connect  TS  to  desired  capacitor.  In 
the  patching  shown  In  Figure  2.2  RT  and  TS  are  determined  from  the  key- 
board. Reset  and  hold  may  also  be  determined  at  the  logic  panel  by 
patching  Dl4  to  R,  015  to  H where  Dlb  and  D 1 5 are  set  by  clock  1 via  the 
logic  board  for  repetitive  operation. 

It  should  be  noted  that  the  output  of  an  operational  amplifier  Is 
the  Inverse  of  the  Input,  l.e. 

Eo  « -AE  jn 
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Figure  2.2A  Logic  Board  Comcor  5000  for  High 
Gain  and  Network  Borrow  Patching 
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where  E0  Is  the  output  voltage,  Ejn  Is  the  Input  voltage  and  A 1$  the  gain 
across  the  summer/integrator. 


I > 


Additional  Inputs  If  needed  may  be  obtained  by  connecting  spare 
resistors  located  on  patchboard  to  SJ  of  the  operational  amplifier. 

Each  of  the  Integrators  has  an  Initial  condition  (1C)  Input  from 
which  IC's  may  be  set  via  potentiometers  or  DAC's.  There  Is  a sign  change 
between  1C  Input  and  the  output  of  the  Integrator. 

Several  static  test  outputs  are  available  via  the  analog  patchboard. 
With  the  analog  In  the  static  test  mode  +100  volts  is  applied  via  these 
terminals  from  which  amplifiers  and  potentiometers  may  be  checked  stati- 
cally. 


Even  numbered  amplifiers  (AOOO  through  A03k)  may  be  patched  as 
high  gain  amplifiers  by  grounding  the  high  gain  connection  located  at  the 
top  of  the  logic  patchboard  thus  eliminating  the  feedback  network.  Input 
networks  may  be  borrowed  from  the  even-numbered  amplifiers  by  patchtng 
the  network  borrow  connections  at  the  top  of  the  logic  patchboard  to 
ground.  See  Figure  2.2A  for  logic  patchboard. 

Each  amplifier  on  the  Comcor  500  may  be  connected  as  a high  gain 
amplifier  whereas  network  borrows  are  available  oniy  on  the  even-numbered 
ampl (flers. 

a 

Multlpl ler 


( 


\ i 
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The  Internal  circuit  and  patchboard  layouts  (stgnel  and  logic)  are 
shown  In  Figure  2.3.  This  unit  performs  division  or  multiplication 
according  to  the  selected  option. 


To  perform  the  division  operation  connect  dividend  analog  Input  to 
Z and  divisor  analog  Input  to  Y.  For  division  tt  Is  also  necessary  to 
connect  Indicated  multiplier  patch  point  on  the  logic  board  to  ground. 


To  multiply,  plus  and  minus  Y Inputs  and  the  +X  Input  must  be 
generated.  It  should  be  noted  for  the  multiplier  shown  minus  X Is 
generated  Internally.  However,  on  some  of  the  multipliers  which  can 
easlty  be  recognized,  mtnus  X is  not  generated  Internally. 

Spare  amplifiers  which  may  be  used  as  Invertors  (M100  to  H123)  are 
available  when  multipliers  are  not  being  used. 


( 


{ 
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and  Outputs 


Figure  2.6  Resolver  2R  Mode  Inputs  and  Outputs 


irWCosi  re  Generator  Inputs  and  Outputs 


Resolver 


Six  operating  options  are  available  when  using 
options,  and  the  necessary  patching,  are  described  as 


the  r*.so1v 
fol lows: 


The??, 


(1)  Sine  and  cosine  mode  - connect  Input  which  must  be  scaled 
1.8  deg/volt  to  0 Input  -0,  +s!n0,  and  +cos0  are  availaf 
outputs.  No  patching  Is  required  on  the  logic  patchboard 
(see  Figure  2.4). 

(2)  1R  mode  - It  Is  necessary  to  patch  1R  on  logic  board  to  logic 
ground.  Shown  In  Figure  2.5  are  necessary  inputs  and  available 
outputs. 

(3)  2R  mode  ~ It  Is  necessary  to  patch  2R  on  logic  board  to 
logical  ground.  Shown  In  Figure  2.6  are  required  Inputs  and 
available  outputs. 

(4)  ROT  Mode  - Rotation  or  transformation  of  axis  mode  requires 
that  ROT  on  logic  board  be  patched  to  logic  ground. 


Inputs  Outputs 

0 Y?  ■ Y.  cos0-X,  sln6 

X.  to  Rj  Input  X,  - Xj  COS0+Y,  slnG 

Yj  to  R2  Input  21  ' 

(5)  Rate  Mode  (continuous  resolution)  - All  of  the  forward  resolu- 
tion modes  (sine,  cosine,  IP.,  2R,  and  ROT)  can  be  extended  to 
Include  continuous  resolution.  An  external  Integrator  Is 
necessary.  To  use  the  Rate  Mode  the  following  Is  necessary: 


OAI  (Digital  to  Analog  switch)  and  RLA  on  logic  board  must  be 
connected  and  RST  must  be  connected  to  computer  RESET  so  that 
the  resolver  starts  In  the  same  quadrant  every  time. 

For  analog  patching  connect  k to  1 Input  and  S001  to  an 
external  Integrator  through  a gain  of  2 as  shown  In  Figure  2.7. 


(6) 


Inverse  Mode 
INV  on  logic 

Inputs 

Xj  Into 

Yj  Into 


(rectangular  to  polar)  - It  is  necessary  to  patch 
board  to  logic  ground. 

Outputs 

R,  Mtan"1  Xl/Yj)  H048 

R2  -R-(A^Y2)  H049 
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Comparator 


The  comparator  compares  two  Input  voltages  and  produces  a logic 
'true'  or  ’false'  depending  on  whether  the  first  Input  Is  greater  or  less 
than  the  second  (or  reference)  Input.  The  logical  ‘not*  of  this  output  Is 
also  aval  able. 

Analog  and  logtc  patchboard  for  a comparator  Is  shown  In  Figure  2,8. 
The  comparator  output  usually  drives  a relay,  as  shown  In  Figure  2.8. 

The  voltage  to  be  compared  Is  applied  to  COO  while  the  negative  of 
the  reference  level  Is  applied  to  B.  When  the  compare  conditions  are 
satisfied  CO  (CO)  on  the  logic  board  change  state.  Either  CO  or  CO  may  be 
used  to  drive  KO  which  switches  the  relay  on  the  analog  patchboard.  An 
external  signal  other  than  CO  or  CO  may  be  used  to  drive  KO  If  desired. 

For  the  relay  outputs,  the  NC  (or  'normally  closed')  position  Is  that  for 
a no  logic  signal  input. 

DA  Switch 

The  OA  switch  provides  a means  of  closing  or  opening  an  analog  signal 
path  with  a logical  variable.  The  analog  and  logic  patchboard  connections 
for  a typical  DA  switch  Is  shown  In  Figure  2.9.  These  are  available  on  the 
Cl -5000  only. 

The  signal  to  be  passed  Is  connected  to  DADO  (analog).'  The  output 
of  the  DA  switch  Is  connected  from  SJ  of  DA  switch  to  SJ  of  an  amplifier. 
Switching  Is  controlled  totally  by  DAO  on  the  logic  patchboard.  Logic 
ground  ('false')  closes  the  switch  while  -6V  opens  the  switch. 

Potentiometer 


Shown  In  Figure  2.10  are  the  schematics  for  the  manual  and  servo  set 
potentiometers.  The  computer  should  ue  In  the  Pot set  mode  when  setting 
each  of  the  potentiometers.  The  servo  potentiometer  Is  set  via  the  key- 
board servo  button  and  DVM  (digital  voltmeter). 

(1)  Place  computer  In  POTSET  mode. 

(2)  Select  pot  to  be  set  via  keyboard.  Pot  number  Is  displayed. 

(3)  Set  DVH  to  value  desired  by  pressing  * on  keyboard  followed 
by  4 digit  potsetting.  Value  displayed  on  OVH. 

(4)  Press  servo  button  and  pot  will  set. 


per£Nr/o  m l rea. 


Figure  2.10  Potent lofnet.r  Patching,  Servo  and  Manual 


Track(i:)-No  Track(f)  Amplifier  . 


T-T  amplifiers  operate  on  the  Comcor  C 1 -5000  as  gain  switches 
(switching  Input  resistors).  The  switch  Is  a DA  switch  controlled  from 
the  logic  board;  l.e.  ground  (false)  DS00  on  logjc  board  results  In  a 
gain  of  10  (T  Input),  amplifiers  100-119  are  TT  amplifiers  while  DS00- 
DS0J9  on  logic  board  are  the  logic  controls. 

Limiter 


The  limiter,  a common  non-llnear  element  In  a system,  Is  available 
directly  on  the  analog  patchboard.  Connect  G of  limiter  to  G of  the  ampli- 
fier to  be  limited  and  LOO  of  Umlter  to  the  output  of  the  amplifier.  With 
the  computer  In  the'POTSET  mode  set  limits  via  the  limit  switches  arid  pots 
located  directly  above  the  analog  patchboard.  With  the  switch  left  set 
lower  limit  with  pot  on  left.  With  the  switch  right  set  upper  limit  with 
pot  on  right.  Monitor  output  of  amplifier  when  limits  are  being  set. 

Function  Switches  (Analog) 

Function  switches  located  to  the  left  of  the  analog  patchboard  are 
double-pole  triple-throw  switches.  Three  different  signals  can  be  patched 
Into  the  L,  C,  and  R positions  and  the  output  at  the  switch  arm  can  be 
selected  by  setting  the  function  switch  to  the  desired  position.  On  some 
switches  A and  B switch  arm  positions  are  available  permitting  simultaneous 
control  of  two  sets  of  function  with  one  switch. 

Trunks 


Trunk  lines  are  available  between  the  Comcor  5000  and  the  Comcor  500 
whlcli  makes  It  possible  to  transmit  signals  between  the  two  machines.  In 
addition.  It  Is  possible  to  slave  the  C I -500  mode  to  the  Ct-5000  mode. 

Trunk  lines  are  also  available  between  the  Comcor  5000  and  the  Angular 
Motion  and  Target  Simulators. 

In  some  cases  ft  Is  desirable  to  use  logic  signals  on  the  CI-5000 
analog  patchboard  (or  vice  versa).  For  this,  trunk  lines  ( I P00- 1 PI  5)  are 
available  between  the  logic  and  analog  patchboards  of  the  Ct-5000. 

2,2.2  Logic  Components 

Single  Input  and  4-lnput  NANO  gates,  delay  flip-flops,  RS  flip-flops 
and  digital  switches  are  available  for  generating  control  logic. 

Delay  time  for  the  delay  flip-flop  Is  determined  by  a rotary  switch 
to  the  right  of  the  analog  patchboard.  Delay  times  available  are  100  usee, 


1 MS,  tO  MS,  100  MS,  t second  and  10  seconds.  Each  switch  also  has  a 
Vernier  scale  which  varies  the  delay  between  the  given  switch  setting  and 
the  next,  thus  covering  the  range  of  0 to  10  seconds  continuously.  To 
activate  the  delay  flip-flop  ground  ENB  and  switch  clock  Input  from  -6V 
to  ground  as  shown  In  Figure  2.11. 
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CLK  I 

_o-q_! 

ENB  | 

0—0 
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DFn 
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0-0 

DFn 

0-0 


Figure  2.11  Delay  Flip-Flop  Logic 


Patchboard 


The  general  purpose  flip-flop  is  shown  In  Figure  2.12  RO  and  SO 
Inputs  provide  direct  flip-flop  control.  These  Inputs  override  any  other 
flip-flop  Input  and  provide  asynchronous  control . A logical  0 (gnd) 
applied  to  either  RO  or  SO  causes  the  FO  or  FO  output  respectively  to  be 
a logical  1 (-6V).  If  a logical  0 Is  applied  to  both  RO  and  SO,  both  FO 
and  FiQ  will  be  logical  1. 


» 

t 


( 


( 


CLK  RO  SO 

O O 0 

GO  FO 

o-  o — o 

Hr  re 

o o — 

Figure  2.12  General  Purpose  Fltp-Flop  Logic 
Patchboard 

The  CLK  terminal  Is  used  for  a clocked  or  asynchronous  signal.  The 
CLK  driver  Is  connected  Internally  to  the  GO,  GO  Inputs.  The  transition 
from  a logical  l to  0 causes  the  flip-flop  to  be  placed  In  a state  given 
by  the  following  J-K  flip-flop  truth  table,  where  J , are  the  present 
states  of  G,  G.  Qn  Is  the  present  output  state  and  Qn+  Is  the  output 
state  after  the  next  clock  pulse. 


i 


2-21 


k“ 

qH+1 

0 

0 

0 

' 1 

0 

I 

0 ; 

1 

1 

— 

■ 

1 ; 

<5* 

Figure  2.13  RS  Flip-Flop  In  J-K  Mode 

Digital  switches  are  different  between  the  two  analog  machines.  Cl- 
5000  switches  are  double-pole,  double-throw  for  switching  logic  signals. 
CI-500  switches  have  single  outputs  which  are  ’true',  'false'  and  pulse 
'true'  for  up,  center  and  momentary  down  positions. 

Control,  sense  and  Interrupt  lines  are  available  between  the 
digital  and  analog  computers.  Control  lines  are  logic  (or  binary)  outputs 
from  digital  to  analog  and  sense  lines  are  logic  Inputs  to  the  digital 
from  the  analog  computers.  Interrupt  lines  are  logic  Inputs  to  the  digital 
which  cause  a program  Interruption  when  the  digital  computer  is  programmed 
to  accept  them.  Interrupt  Inputs  should  be  clocked  pulses;  cqntrol  and 
sense  lines  are  logic  levels. 

DAC  and  ADC  communication  between  digital  and  analog  computers  Is 
available,  Patching  for  Input/output  is  placed  on  the  analog  patchboard, 
but  S/H  control  for  ADC's  and  DAC  update  control  Is  performed  by  patching 
special  control  lines  on  the  logic  board.  Each  special  control  line 
controls  2 ADC's  or  2 DAC's  (4  each  on  the  Cl-500)  and  Is  operated  by  a 
pulse  on  the  line.  All  DAC's  are  of  the  multiplying  type  whereby  an 
analog  Input  signal  at  the  DAC  output  causes  the  digital  and  analog  signals 
to  be  multiplied  together.  If  required,  the  analog  signal  may  be  constant 
at  +100  volts  to  give  a unity  multiplier. 

External  Integrator  control  and/or  repetitive  runs  are  available 
via  the  clock  1 Input  located  on  the  logic  board.  Clock  1 Is  driven  by 
clock  frequencies  shown  on  the  logic  patchboard.  Counters,  which  count 
number  of  clock  pulses  while  the  computer  is  In  the  reset,  hold,  and 
compute  modes  are  located  on  the  console  in  respective  order.  For  Clock 
1 to  be  active  the  computer  must  be  in  the  compute  mode  with  the  clock 
button  set.  The  associated  run,  hold  control  signals  (Rl , HI)  are 
generated  via  the  logic  patchboard  and  may  be  used  to  control  Integrator 
states.  See  Figure  2.14  for  setup. 
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Figure  2.14  Logic  Board  for  External  Clock  Control 


Comcor  500 


Analog  and  logic  patchboards  of  the  C I —500  computer  contain  connections 
which  are  not  Indicated  on  some  patchboards.  These  connection  points  are 
Indicated  as  follows.  It  should  be  noted  that  C I —500  patchboard  component 
numbering  Is  in  octal  format: 

DAC's  and  ADC's  appear  on  the  analog  board  as  shown  In  Ftgure  2.15. 


In 


i 


A A 

202  G 203 


i 


DAC  0 DAC  1 


DAC  2 DAC  3 


DAC'S 


In  Out 


DAC  31 


AAA 


i | ■ Si 

ADC0  A0C1  ADC2  ' ADC 3 ADC31 


ADC'S 


Figure  2.15  C I -500  DAC  and  ADC  Patch  Connections 


On  the  C 1 -500  logic  board,  the  control  Inputs  for  ADC  Sample/Hold 
and  DAC  update  are  Indicated  In  Table  III.  There  are  no  special  control 
lines  for  these  functions.  Each  control  Input  operates  on  a group  of 
four  corresponding  ADC's  or  DAC's. 

Patching  for  the  two  C I -500  Interface  clocks  Is  given  In  Table  IV. 

For  these  clocks,  Input  control  points  are  available  at  the  analog  console, 
as  are  status  monitoring  points.  Clock  0 on  the  Cl-5000  Interface  has 
only  pulse  out  available  on  the  logic  patchboard,  plus  an  Input  for  an 
external  clocking  rate  to  replace  the  Internal  rate  of  1 microsecond. 


r 


Logic 

Address 

Clock  1 
Function 

Logic 

Address 

Clock  2 
Function 

T 020 

RESET 

IN 

* 

RESET 

IN 

T 021 

STOP 

(N 

* 

STOP 

IN 

T 022 

RUN 

IN 

• * 

RUN 

IN 

T 023 

RESET 

OUT 

* 

RESET 

OUT 

T 02k 

STOP 

OUT 

* 

STOP 

OUT 

T 025 

RUN 

OUT 

i 

RUN 

OUT 

T 026 

PULSE 

OUT 

* 

PULSE 

OUT 

TABLE  IV  Cl -500  Clock  Patching 


♦These  patchpofnts  are  unmarked,  but  are  located  In  the 
white  painted  area  occupying  the  six  upper  leftmost  pairs 
of  holes  (l.e.  In  2 columns  and  three  rows). 
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Multipliers  on  the  Comcor  500  require  external  amplifiers.  The  Input 
amplifier  Is  a sign  change  only  while  the  output  amplifier  Is  high  gain. 

The  high  gain  amplifier  Is  generated  by  patching  HG  on  the  logic  patch- 
board for  the  amplifier  used. 

To  obtain  Integrator  mode  control  on  the  Comcor  500  the  logic  board 
must  be  patched  as  shown  In  Figure  2.15A. 


8 6 

HG  .Si  TS  .9 


R H 


Figure  2.15A  Integrator  Mode  Control  Logic 


Other  operation  $ on  the  Comcor  500  analog  and  logic  patchboards  are 
similar  to  the  Comcor  5000  operation. 

For  non-linear  functions,  diode  function  generators  are  available 
on  the  Comcor  5000  (servo-set  - 10  segment)  and  vh  Com -or  500  (hand-set 
- 10  segment).  For  a detailed  description  of  operation  the  reader  Is 
referred  to  the  CI-225A  Manual. 

2.2.3  Summary 

Shown  below  are  tables  summarizing  component1?  of  the  Comcor  500  and 
Comcor  500.  Also  given  Is  trunkline  identification  between  the  Comcor  5000 
and  Comcor  500. 


ANALOG 


5000 

500 

Amplifiers 

120 

40  (0-4?) 

Integrators 

48 

20 

Multfpl iers 

48 

24 

Resolvers 

8 

4 

Comparators 

20 

8 

Relays 

16 

8 

DA  Switch 

16 

0 

Potentiometer 

Handset 

24 

24 

Servo-set 

100 

48 

T-T  Amplifier 

20 

0 

Limiters 

16 

0 

Function  Switch 

20 

8 

ADC's 

32 

32 

DAC's 

Trunks 

500—5000  (31) 

32 

32 

500 — Table  Room  (10) 

(T140-T141) 

5000— Table  Room(27) 

(T066-T093) 

LOGIC 

* 

5000 

500 

' 

' 1 

Hand  Gates 

1 -Input 

24 

18 

2-Input 

24 

16 

4- Input 
Ft  Ip-Flops 

24 

18 

Oelay 

8 

3 

R-S 

32 

16 

♦Interrupts 

16 

16 

Control  Lines 

16 

16 

Sense  Lines 

16 

16 

Digital  Switch 

10 

8 

♦The  XDS  Slgroa  5 has  only  16  interrupts  each  of  whtch  nay  be 
triggered  from  either  the  Coacor  5000  or  Concor  500. 
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5000 

500 

T 032 

T 040 

T 033 

T 04! 

T 034 

T 042 

T 035 

T 043 

T 036 

T 044 

T 037 

T 045 

T 038 

T 046 

T 039 

T 047 

T 040 

T 050 

T 041 

T 051 

T 042 

T 052 

T 043 

T 053 

T 044 

T 054 

T 045 

T 055 

T 046 

T 056 

T 047 

T 057 

T 048 

T 060 

T 049 

T 061 

t 050 

T 062 

Y 05 1 

T 063 

T 052 

T 064 

7 053 

T 065 

T 054 

T 066 

T 055 

T 067 

T 056 

T 070 

T 057 

T 071 

T 058 

T 072 

T 059 

T 073 

T 060 

T 074 

T 061 

T 075 

T 062 

T 076 

T 063 

T 077 

TABLE  V Trunkline  Identification  between 
Coaacor  5000  and  Comcor  500 
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The  following  Is  a list  of  Comcor  5000  power  supply  addresses  and 
values: 


P500 

+100 

P508 

+20 

P501 

-100 

P509 

-20 

P5Q2 

+150 

P510 

-20 

P503 

-150 

P511 

+18 

P504 

+20 

P512 

-18 

P505 

-20 

P513 

-6 

P506 

+150 

P519 

ROAC 

P507 

-150 

The  following  Is  a list  of  Comcor  500  power  supply  addresses  and 
values: 


P500 

+100 

P516 

-6 

P501 

-100 

P517 

-18 

P502 

+150 

P520 

+18 

P503 

-150 

P521 

-6 

P50J* 

+20 

P522 

+26 

P505 

-20 

P523 

-2k 

P506 

+20 

P524 

-2k 

P507 

-20 

P525 

RDAC 

P510 

+150 

P526 

-18 

P5H 

-150 

P527 

+18 

P512 

+20 

P530 

-6 

P5I3 

-20 

PS1A 

+18 

P515 

-18 
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As  Illustrated  In  Figure  2.1,  the  digital  computer  Is  Interfaced  to 
Astredata  Cc„icor  Cl -500  and  C I —5000  analog  computers.  Data  Is  transferred 
through  the  Interfaces  via  DAC's,  ADC's  and  control  and  sense  (i.e. 
discrete  or  binary)  lines.  Each  interface  has  one  or  more  real-time  clocks 
which  are  accessible  at  both  sides  of  the  Interface  to  permit  time  synchroni- 
nlzation  of  analog  and  digital  elements.  Inputs  to  the  Sigma  5 external 
Interrupt  system  are  also  available  at  both  analog  consoles.  Table  VI 
lists  the  data  transfer  elements  In  each  Interface.  Note  that  the  special 
control  lines  In  the  CI-5000  Interface  are  for  use  In  controlling  simul- 
taneous DAC  updating  or  ADC  sampling. 

2.3.1  Signal  Transfer 

Both  Interfaces  are  connected  to  the  Sigma  5 DIO  channel  and  all  data 
and  logic  signals  may  be  transmitted  by  appropriate  ^rlte  Direct'  or  'Read 
Direct1  commands  (see  Reference  1),  However,  DAC's  and  ADC's  In  the  Cl -500 
Interface  may  also  transmit  data,  under  external  Interrupt  control,  to  and 
from  memory  directly  by  means  of  the  high  speed  DMA  channel.  Continuous 
operation  In  the  DMA  mode  may  be  obtained  through  use  of  the  8 pre-set 
registers  of  the  DMA  In  combination  with  the  'chain'  option. 

Discrete  signals  may  be  transmitted  singly  or  In  groups.  Reference  2 
contains  details  of  the  various  commands  and  Interface  1/0  modes. 

2.3.2  DAC  and  ADC 

Digital-analog  and  analog-digital  conversion  uses  15  bit  wide  digital 
data  which  represents  an  analog  range  of  +100  volts.  Both  DAC's  and  ADC's 
have  optional  fixed  and  floating  point  modes;  In  floating  point  mode  0AC 
inputs  and  AOC  outputs  are  In  standard  floating  point  format  In  the  range 
+1.0  and  In  this  case  the  data  occupies  an  extra  8 bits  to  allow  space  for 
the  floating  point  exponent.  In  fixed  point  mode  the  digital  value  from 
0AC  or  to  AOC  occupies  a 32  bit  word  with  the  most-significant  15  btts 
representing  a binary  fraction  In  the  range  -1  to  +(l-2"'  v with  negative 
numbers  In  2-complement  form. 

DAC's  are  of  the  multiplying  type  wherein  the  analog  output  Is  formed 
by  multiplying  the  converted  digital  signal  by  an  analog  signal  which  Is 
Input  at  the  analog  patchboard.  A unity  analog  multiplier  may  be  obtained 
by  patching  a +100  volt  signal. 

For  ADC's  which  terminate  on  the  CI-5000  (addresses  0-3!) , a sign 
Inversion  occurs  on  data  transfer.  Positive  voltages  Input  to  the  ADC 
are  read  as  negative  data  In  the  digital  computer.  On  the  Cl -500  the 
DAC's  (addresses  32-63)  Invert  so  that  a positive  value  output  by  the 
digital  computer  produces  a negative  voltage  If  the  analog  multiplier  Is 
positive.  AOC  conversion  on  the  Cl-500  ^addresses  32-63)  and  0AC  con- 
version on  the  CI-5000  (addresses  0-31)  do  not  Involve  a sign  Inversion. 

<: 
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2. 3 <■  3 Control  and  Sense  Lines 


Control  lines  may  be  set,  reset  or  pulsed  singly  and  may  be  set  or 
reset  in  a group  of  32  or  16  (for  C ? —5000  and  C i -500  respectively).  The 
status  of  the  32  group  1 (C I -5000)  control  lines  may  also  be  read  as  a 
group.  Sense  lines  may  be  read  singly  or  In  the  same  groups  of  32  and 
16. 


Special  control  lines  are  available  (on  the  C I -5000  only)  for  con- 
trol of  DAC  transfer  and  ADC  sample/hold.  If  the  control  inputs  to  DAC 
and  ADC  are  patched  then  DAC  conversions  are  only  updated  when  a pulse 
Is  sent  to  the  control  and  ADC  sample  or  hold  mode  is  dependent  on  the 
control  input  being  reset  or  set.  If  the  control  inputs  have  nothing 
patched  to  them,  then  DAC  updating  and  ADC  S/H  occurs  whenever  an  1/0 
operation  is  requested  by  the  digital  machine.  On  the  Cl -5000  there  are 
16  DAC  and  ADC  control  Inputs,  each  operating  on  a pair  of  DAC’s  or  ADC's. 
The  C I “500  has  8 DAC  and  8 ADC  control  inputs,  so  that  each  input  controls 
k DAC's  or  A ADC's;  the  C i —500  DAC  and  ADC  control  Inputs  are  not  access- 
ible dlractly  through  the  A/D  interface,  but  only  through  logic  signals 
patched  at  the  logic  patchboard. 

2.3.^  Interrupts 

The  Sigma  5 digital  computer  contains  16  external  Interrupts,  numbered 
0 to  15  external  group  2.  interrupt  0 Is  permanently  connected  to  the 
DVM's  (Digital  Voltmeters)  of  the  two  analog  computers  and  Is-not  available 
at  the  logic  patchboards.  External  Interrupts  1 to  15  are  available  at 
both  analog  consoles  and  are  paralleled  into  the  Sigma  5 but  Interrupt  1 Is 
wired  to  the  DMA  channel  and  should  not  be  used  by  a foreground  program. 
Interrupt  response  processing  in  the  digital  machine  has  no  method  of 
Identifying  the  source  of  an  external  Interrupt,  and  thus  care  should  be 
exercised  to  avoid  confusion  in  the  use  of  the  external  Interrupts.  In 
addition,  the  graphics  display  unit  Interface  (7611)  Is  wired  Into 
Interrupts  6 and  7 and  these  Interrupts  should  not  be  used  at  the  analog 
consoles  If  the  graphics  terminal  is  in  use  by  any  (foreground  or  background) 
program.  Furthermore,  the  HRBM  monitor  Control  Task  uses  external 
Interrupt  15  and  thus,  this  Interrupt  should  not  be  used. 

2.3.5  Clocks  (Interval  Timers) 

The  Cl -5000  Interface  contains  an  Interval  timer  (Clock  0)  with  the 
following  characteristics: 

(a)  Selectable  Internal  and  external  modes 

(b)  Single  or  repeating  countdown  (internal  mode) 
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(c)  Internal  count  rate  of  I megahertz  (1  microsecond). 

(d)  Maximum  Internal  cycle  count  65 s 537 * 

(e)  External  mode  Input  at  Cl -5000  logic  patchboard. 

(f)  Control  from  digital  computer  DIO  channel  only. 

The  C I —500  Interface  contains  two  Interval  timers  (Clocks  1 and  2) 
with  the  following  characteristics: 

(a)  Internal  mode  only. 

(b)  Repeating  countdown  only. 

(c)  Courft  rate  100  kilohertz  (10  microseconds). 

(d)  Maximum  cycle  count  65,537* 

(e)  Control  from  digital  computer  DIO  channel  or  C I -500  logic 
patchboard. 

Control  operations  consist  of  resetting,  setting,  Initializing, 
running,  stopping  and  reading  the  clocks.  In  addition,  the  status  of 
Clocks  1 and  2 may  be  determined  at  the  Cl -500  logic  patchboard. 

* 

2.4  Angular  Motion  Simulator 

2.4.1  Hardware  Characteristics  and  Axes  Systems 

The  angular  motion  simulator  consists  of  a hydraulic,  servo-driven 
platform  mounting  with  three  gimballed  rotational  degrees  of  freedom 
about  mutually  perpendicular  axes  and  associated  control  and  hydraulic 
systems.-  Dynamic  Inputs  to  the  control  system,  originating  from  a real- 
time dynamic  simulation,  enable  hardware  components  mounted  on  the  plat- 
form to  experience  real  In-service  rotational  rates  and  posttlons. 

Table  VII  shows  nominal  performance  specifications  of  the  motion  simulator. 

GImbal  order  Is  roll,  yaw,  pitch  with  roll  the  Inner  axls>  Thus, 
a convenient  Euler  angle  representation  for  defining  orientation  of  body 
axes  of  the  hardware  components  relative  to  an  Inertial  reference  Is  the 
non-orthogonal  ordered  set  4',  0 where  the  angles  are  conventionally 
def 1 ned 

Inputs  to  the  simulator  consist  of  three  volt«ages  which  determine 
the  angular  positions  of  the  three  axes.  Plckoff  output  voltages  are 
available  which  yield  the  actual  position  of  each  axis  at  any  Instant. 
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TABLE  VII 

Angular  Motion  Simulator 
Performance  Specifications 

with  5.0  lb-ln-sec2  load 
0 930  lb-ln2) 


Roll 

(red) 

Yaw 

(yellow) 

Pitch 

(blue) 

Max.  Acceleration,  exceeds 

50,000°/sec2 

25,000°/sec2 

12,000°/sec2 

Max.  Velocity,  exceeds 

700°/sec 

400°/sec 

200°/sec 

Min.  Velocity,  less  than 

0.0004°/sec 

O.OOOWsec 

0.0004°/sec 

Displacement 

+120° 

+120° 

+120° 

Position  Accuracy  (Maximum) 

+0.053° 

+0.053° 

+0.053° 

Repeatability,  less  than 

+0.005° 

+0.005°  * 

+0.005° 

Readout  accuracy 

+0.12° 

+0.12° 

+0.12° 

Command  $ 1 gna 1 Sea  1 1 ng 
(typical) 

1.80°/volt 

1.80°/volt 

l.80°/volt 

Frequency  Response 

27.5  cps 

17.6  cps 

12.8  cps 

Nominal  scalings  between  angular  position  and  voltage  are  shown  In  Table 
VII.  These  may  easily  be  changed  to  produce  Increased  sensitivity  at  the 
cost  of  reduced  displacement  ranges  by  varying  the  amplifier  Input 
resistance.  Maximum  Input  voltages  are  +100  volts. 

To  define  the  sign  conventions  of  the  glmbal  angles  <{>,  0 (roll, 

yaw,  pitch)  consider  the  set  of  Cartesian  axes  X,  Y,  Z with  origin  at 
the  center  of  gimballed  rotation.  The  X axis  Is  horizontal,  positive  de- 
fined geographically  as  pointing  to  the  south  wall  of  the  laboratory,  the 
Y axis  Is  positive  vertically  upwards  and  the  Z axis  to  the  right  when 
facing  In  the  +X  direction.  Positive  voltages  Into  the  table  control 
system  produce  glmbal  angles  corresponding  to  an  ordered  set  of  Euler 
angles  In  the  above  order  {$,  4>,  0)  which  locate  a set  of  axes,  fixed  In 
the  Inner  glmbal,  relative  to  X,  Y,  Z. 

A common  set  of  axes  used  In  aerospace  simulations  Is  that  obtained 
by  rotating  the  X,  Y,  Z set  n/2  about  X (I.e.  Z vertically  downwards). 

In  this  case  the  tp  Euler  angle  Input  voltage  to  the  simulator  requires  to 
be  Inverted  In  sign  for  correct  response  of  the  yaw  glmbal.  0 and  are 
unchanged. 


2. 4. 2 Phase  Compensation 

Figures  2.16  and  2.17  Illustrate  the  frequency  response  of  the  pitch 
and  yaw  channels  of  the  motion  simulator.  Occasionally  the  phase  lag 
exhibited  is  not  acceptable  and  the  application  of  a lead  filter  to  the 
Input  signals  Is  necessary.  A network  with  a frequency  domain  response 
given  as: 

f (s)  — f t X s + 1 | 2 

where  tj  ■ .0222  and  12  * .0071,  has  been  found  to  give  adequate  compen- 
sation to  signals  which  contain  frequencies  to  about  4-5  Hz.  The  response 
of  the  above  filter  Is  Included  In  Figures  2.16  and  2.17  and  an  analog 
computer  patching  diagram  which  Implements  the  filter  Is  shown  In  Figure 
2.18.  Since  this  filter  Is  essentially  a differentiating  network,  Input 
signals  should  not  be  derived  from  OAC's  without  some  Intervening  smoothing 
to  give  a continuous  Input  signal  with  continuous  first  derivative. 

2.4.3  Simulation  Hardware  Protection 

When  running  hardware-ln-the-loop  simulations  which  Include  gyro- 
mounted  seekers  or  similar  units  on  the  Angular  Motion  Simulator,  It  Is 
usually  wise  to  protect  the  hardware  components  by  an  automatic  run-abort 
mechanism  based  on  gyro  precision  angles  approaching  a maximum  value. 

This  Is  achieved  by  feeding  the  gyro  angles,  which  are  usually  available 
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Figure  2.18  Analog  Computer  Mechanization  of  Lead  Filter  I 
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at  test  points,  to  one  of  the  analog  computers,  forming  the  total  angle 
and  comparing  this  with  a pre-determlned  limit.  If  the  limit  Is  reached, 
a relay  Is  used  to  set  a logic  true  pulse  which  generates  an  Interrupt 
signal  to  the  digital  control  program  enabling  the  simulation  to  be 
aborted  before  the  gyro  reaches  Its  mechanical  stops.  The  analog  mechani- 
zation of  this  circuit  Is  shown  In  Figure  2.13,  assuming  gyro  angle 
components  0S,  $s. 

2.5  Target  Simulators 

The  target  simulators  consist  of  laser  and  optical  beam  generators 
which  project  either  a laser  spot  or  slide  view  of  a target  respectively 
onto  the  rear  of  a transulucent  screen.  The  ground  glass  translucent 
screen  causes  a diffused  Image  to  be  visible  to  hardware  mounted  on  the 
Angular  Motion  Simulator.  The  laser  simulates  the  reflection  of  a laser 
designator  fronf  a target,  and  the  optical  contrast  generator  simulates 
the  view  of  a target  as  seen  by  an  on-board  video  receiver.  Position  of 
the  simulated  target  on  the  screen  Is  controlled  by  deflecting  the  beam 
by  means  of  one  (optical  simulator)  or  two  (laser  simulator)  mirrors 
which  rotate  about  axes  In  the  plane  of  the  mirror  with  two  degrees  of 
freedom.  The  position  of  the  simulated  target  on  the  screen  can  be  con- 
trolled to  simulate  the  actual  line  of  sight  angles  between  missile  and 
target  when  the  missile  seeker  Is  mounted  on  the  angular  motion  simulator. 
Target-missile  closure  Is  simulated  by  Increasing  the  spot  size  In  the 
laser  simulator  and  zooming  the  slide  display  In  the  optical  simulator. 

The  simulator  configuration  In  the  laboratory  Is  shown  dlagrafhmatlcally 
In  Figure  2.20.  Voltages  to  control  mirror  deflections,  laser  beam 
aperture  position  and  zoom  lens  position  may  be  set  manually,  derived 
from  computer  outputs  (analog)  or  other  external  sources. 

2.5.1  Angular  Relationships 

Displacement  of  the  simulated  target  display  to  represent  the  rela- 
tive motion  of  missile  and  target  Is  governed  only  by  the  relative  dis- 
placement coordinates  between  the  missile  and  the  target.  The  mathematics 
underlying  the  derivation  of  mirror  deflection  angles  are  given  In 
reference  3,  but  will  be  summarized  hera. 

Axes  systems  used  In  the  derivation  are  as  follows: 

(a)  The  Cartesian  Inertial  system  translating  (or,  Instantaneously 
coincident)  with  the  missile,  with  the  CG  as  origin.  Missile 
angular  orientation  Is  defined  relative  to  these  axes. 

(b)  The  Cartesian  axis  system  defined  for  the  Angular  Motion  Simu- 
lator (Z  vertically  downwards  and  Y horizontal). 


(c)  Two-dimensional  axes  y,  z defining  target  display  displacement 
on  the  screen,  as  shown  In  Figure  2.21. 

Since  the  screen  Is  approximately  5 feet  by  5 feet,  the  range  of  line 
of  sight  angles  which  can  be  simulated  are  determined  by  the  requirement 
of  keeping  the  display  on  the  screen  without  reaching  the  limits  of  mirror 
deflection.  It  is  particularly  Important  that  the  mirror  does  not  reach 
Its  limit  of  travel  If  the  seeker  uses  target  line  of  sight  rate.  In 
order  to  select  and  control  the  range  of  line  of  sight  angles  which  are  to 
be  displayed  during  a simulation,  It  is  assumed  that  the  axis  system  (a) 

Is  Inclined  to  that  of  (b)  at  an  angle  which  may  be  expressed  as  a component 
0|  about  Y in  (b)  and  i|>|  about  the  Z In  (a).  Actual  target  missile 
relationships  are  Illustrated  In  figure  2.22,  where  the  line  of  sight 
angular  components  are  functions  of  the  missile-target  separation  distances 
AX,  AY,  AZ  and  are  given  by: 

• 

-AZ 

tan  8lo$  - AX 
AY 

tan  *L0S  VaxW 

From  these  definitions,  the  displacement  of  the  target  display  on  the 
screen  l s (from  Reference  3) : 


y • Lsy  cos  (8,  cosi{> | ) cos  | cosO|  n<;)tan(-j>|  cos9L0S  ♦ 
COS9|  COSlJ»|  cos'to  | COSty  I ♦ 0|_qcJ" 


cos(6|  eosfri)  tan (8 1 cosifri  + 6ms) 

CQS0 | COS^j 

Is  the  distance  frrxa  the  display  screen  to  the  Angular  Motion 
center  of  rotation. 

angles  of  deflection  of  the  laser  or  optical  beam  onto  the 
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where  Is  the  distance  from  mirror  number  1 to  the  screen. 

For  the  particular  case  of  i|>|  ■ 0 In  the  above  equations  for  y and 
z,  we  have 


AY 

AXcos8|  + AZsinQj 


Lr_  AX  tanflj  - AZ 
ST  AX  AZ  tan 6 1 


In  the  case  where  ^ j " 0 we  have  a conmon  configuration  wherein  the 
actual  Inertial  axes  coincident  with  the  missile  are  related  to  the 
laboratory  axes  (whose  X axis  Is  the  horizontal  line  from  Motion  Simulator 
to  screen)  by  a rotation  8.  about  their  common  Y axis.  By  varying  0|,  a 
range  of  6^  is  selected  for  display  on  the  target  simulator  screen. 

For  the  general  case  of  6 , ij>.  values  the  Inputs  to  the  Angular 
Motion  Simulator  0-j-.  are  (from  Reference  3): 


6T  ■ tan”1  f-bjj  costyj  tanflj  - bt?  slmfrj  tan6|  - b??^ 


b|  | cosi|i|  + bjj  slni^i  + bjj  tan0j  f 


*r  “ t®""1  J h l.-ilrii JL..b,i2^ii 


,/\  - <b| | slmfij  ♦ bj2  cos<»|)2  J 


♦t  " tan~'  j ^3iJ»iU..l432 

b2l  Un  I 4 b22 


where  b j j are  the 
missile  body  axes 

b|  I “ COS\J>£ 

b2,  - sln*E 

b31  * cos*E 
b,2  - sln*E 

b22  * cos*E 


elements  of  the  transformation  matrix  from  Inertial  to 
through  Euler  angles  0E,  tE,  +£.  Thus, 

cos0E 

sln0E  - co$$E  sln^-  cos0E 
slnOg  ♦ sln$E  sIn<;E  cos©E 

cos*E 
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t>32  ■ -sln^  cos^E 
b|3  ■ “cos^E  sin0£ 

- stn«j>£  cos0£  + cos4)g  sfn^E  stn0E 
^33  " C0S(^e  cos®£  " sln$E  sln^E  stn0E 


For  the  case  where  if»|  » 0,  the  Angular  Motion  Simulator  angles  take 
the  form: 

*t|  . “9Et6l 


* VE 

I— Q 


*1. 
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2.5.2  Simulation  of  Target  Closure 

Closure  of  the  range  between  target  and  missile  Is  simulated  for  laser 
seekers  by  Increasing  the  reflected  spot  size,  and  hence  the  reflected 
laser  energy,  proportionally  according  to  the  Inverse  of  the  slant  range 
between  target  and  missile.  Laser  spot  size  Is  varied  by  passing  the  laser 
beam  through  a convergent  lens  and  then  through  a moveable  knife-edge 
circular  aperture,  which  produces  a spot  diameter  ratio  of  80:1  (38  db) 
between  Its  closest  and  furthest  positions  from  the  convergent  lens.  A 
further  30  db  of  energy  attentuatlon  Is  available  In  10  db  steps  by  means 
of  filter  attenuators  placed  ahead  of  the  convergent  lens.  Figure  2.23 
shows  the  variation  of  laser  energy  density  at  the  screen  with  aperture 
position  which  can  vary  from  2.25  inches  to  22.0  Inches  from  the  lens. 

Aperture  position  Apr  as  a function  of  slant  range  from  missile  to 
target  Is  given  by  the  following: 

R5R  - /&X*  + AY2  + AZ2 
then  aperture  Input  voltage  Ap^  Is 

*?C  " KAF  when  RSR  > R5Tf  t 
“ KAN  when  ^SR  * SjpLT 
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APC  * MaRSR  + 


where  M„  « -Dd(Ka„  - KftC) 

a -pR — ANTTr^FTi 

rSTRT  ^r  ' > 


DdKau  - K 


r^F 


substituting  for  Ha  and  Ba  gives 

Apc  " & — , §SR_  *KAF  ' KA|P'  ■ £aF  * KAN 

dr  rstrt  dR 

where  K ® 62 *k  Is  the  maximum  voltaci  and  * -9^.0  Is  the  minimum 
voltage.  Is  the  ratio  of  largest  to  smallest  spot  size  taking  the 
smallest  as  unity,  t.e.  «=  80  and  Rctdy  Is  the  slant  range  at  which  the 
spot  size  Is  assumed  to  start  Increasing.  Note  that  the  above  expressions 
do  not  depend  on  the  absolute  values  of  Rs«  and  Reypy  hut  only  their  ratio. 
Also,  the  values  of  K/^p  and  are  given  for  +100  volt  analog  range.  For 
a digital  program  using  the  floating  point  option  on  a DAC  these  values 
are  divided  by  100. 


In  the  case  of  the  optical  contrast  simulator,  target  range  closure 
fs  simulated  by  the  use  of  a high-speed  zoom  lens  with  a magnification 
ratio  of  50:1  and  the  capability  of  simulating  closing  speeds  of  up  to 
1600  ft/sec.  Input  voltage  to  the  zoom  control  unit  varies  from  zero  to 
100  volts  and  is  given  by  the  following  polynomial  fit  to  measured  data: 

Zv  - -25.71^6  + 3!.43604Xm  - S.39579^2  + 0.3^6925XH3 
when  XM  < 7.075 

Zv  * 18.115  + 5.917363Xh  - 0.2292095XM2  + 4.77775*lO'3XH3 
-3.8J2966*!0“5Xm^  when  X„  >,.  7.075 

where  Zy  Is  the  zoom  control  voltage  and 

XH  “ &Z  “ — ^5Z — 

SR  Ax2+ay2+az2 


and  R Is  the  range  at  which  the  magnification  Is  i.0. 
Si 
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The  process  of  zooming  the  optical  contrast  lens  results  In  a 
rotation  of  the  Image  on  the  display  screen.  The  rotation  Is  compensated 
for  by  adding  a correction  term  to  mirror  angles  (both  mirrors,  If  both 
are  being  used)  obtained  by  Interpolation  In  tables  for  the  6^  correction 
foc  and  the  4)mc  correction  f^  as  functions  of  the  zoom  voltage  Zy.  These 
data  are  shown  In  Table  VI 51.  Addition  of  the  corrections  to  both  mirrors 
Is  necessary  because  of  the  opposite  polarity  of  the  mirrors;  the  object 
Is  to  correct  the  zoom  distortion  by  modifying  mirror  I,  but  to  eliminate 
the  corrections  from  the  laser  deflection  by  Inserting  the  Inverse 
corrections  In  mirror  2. 

Note  that  both  optical  and  laser  simulators  will  only  be  used  simul- 
taneously when  using  ‘dual-mode  seekers'.  Target-missile  motion  will 
normally  be  simulated  using  mirror  1,  and  mirror  2,  which  has  a range  of 
motion  of  +3°  In  each  channel,  will  be  used  for  second-order  effects  such 
as  system  noise  and  spot  jitter. 

2.5.3  Hardware  Characteristics 

The  laser  target  simulator  employs  a pulsed  Q-swItched,  neodymlum- 
YAG  laser  operating  at  either  10  or  20  pulses/second.  Physical  layout  of 
the  simulator  Is  shown  In  Figure  2.2k.  The  circular  glrnballed  mirrors 
are  Identical;  each  mirror  has  a range  of  deflection  between  mechanical 
stops  of  +15°  In  both  degrees  of  freedom.  However,  for  the  physical  layout 
shown,  the  laser  or  light  beam  Is  deflected  to  the  edge  of  thp  screen  by  a 
+9°  deflection  of  mirror  number  l and  for  the  laser  beam  to  Impinge  on 
mirror  l,  mirror  number  2 Is  limited  to  deflections  of  +3°*  As  a result, 
main  control  of  the  spot  deflection  Is  exercised  through  mirror  number  1 
and  mirror  number  2 Is  used  to  simulate  laser  spot  positional  noise 
and  jitter  and  to  oppose  corrections  io  mirror  I required  by  the  zoom  con- 
trol of  the  optical  contrast  simulator. 

The  energy  density  of  the^laser  spot  on  the  screen  may  range  from 
1.2*10”^  to  ^.3*10-9  joules/cm‘  . Further  details  of  laser  performance, 
calibration  and  controls  are  given  in  Reference  k. 

The  optical  contrast  target  simulator  consists  of  a photographic 
slide  projector  equipped  with  a high-speed  zoom  lens  which  has  a dynamic 
range  of  50:1  In  Image  magnification.  Magnification  at  the  smallest  zoom 
setting  Is  4X,  giving  about  200X  at  the  maximum  zoom  position.  Image 
size  Is  maintained  at  17  Inches  x 17  Inches  on  the  screen.  Screen  Illum- 
ination Is  greater  than  k ft-cand!es  at  all  zoom  positions  for  an  Input 
power  to  the  projector  lamp  of  about  2 kilowatts  (lamp  rated  at  2.5  kilo- 
watts). The  time  required  to  zoom  the  lens  over  Its  full  ranee  of  travel 
Is  less  than  3.9  seconds  permitting  the  simulation  of  closing  speeds  up  to 
•bout  1 600  ft/second. 
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The  projector  employs  a 2.5  Kilowatt  Xenon  lamp  with  an  emission 
spectrum  covering  wavelengths  from  0.2  microns  to  1.1  microns.  A series 
of  optica!  filters  Is  Inserted  between  the  lamp  and  slide  to  reject  the 
haat  producing  near  infra-red  radiation  at  wavelengths  above  0.72  microns. 
The  heat  is  carried  away  by  cooling  fans.  An  alternative  filter  is 
available  which  permits  radiation  to  0.81  microns  to  pass  through  to  the 
screen.  Care  should  be  exercised  to  limit  the  time  at  high  magnifications 
to  avoid  damage  to  the  slide  when  using  this  alternative  filter.  Neutral 
density  filters  are  available  for  mounting  on  the  zoom  lens  aperture  to 
simulate  various  levels  of  atmospheric  haze.  Further  details  of  set-up, 
use  and  calibration  of  the  optical  contrast  simulator  are  given  In 
Reference  5- 

Voltage  inputs  to  Mirrors  1 and  2,  laser  spot  size  controller  and 
zoom  position  have  been  expressed  as  being  In  the  rangs  +.1 00  volts.  This 
Is  for  analog  computer  compatibility  purposes;  a buffer  amplifier  with  a 
gain  of  0.1  Is  Inserted  In  the  input  circuitry  to  the  hardware  since  these 
components  operate  on  +10  volts.  The  position  of  any  of  the  above-mentioned 
hardware  components  may  be  obtained  at  any  time  by  reading  the  output 
voltage  of  the  desired  unit.  This  voltage  Is  in  the  range  +10  volts. 

Since  the  mirrors  have  a restricted  frequency  response,  as  shown  in  Figures 
2.25  and  2.26,  It  is  often  necessary  to  monitor  the  output  of  Mirror  1 
to  ensure  that  it  does  not  diverge  too  greatly  from  the  input  waveform. 

Too  great  a divergence  may  result  In  unacceptably  large  errors  In  spot 
motion.  Voltage-angle  scale  factors  are  non-linear  for  Mirrors  1 and  2; 
their  variation  Is  shown  in  Figures  2.27  and  2.28.  It  should'be  noted 
that  a positive  voltage  In  the  mirror  pitch  control  axis  (6^)  produces  a 
negative  mirror  angle,  which  Is  a downward  deflection  of  the  screen  Image. 

A positive  yaw  voltage  produces  a positive  yaw  mirror  angle  (4'rnc)  which  Is 
a leftward  movement  of  the  screen  Image  when  viewed  from  the  Angular  Motion 
Simulator. 

2.5.**  Control  and  Patch  Panels 

The  control  console  for  the  target  simulators  contains  the  following: 

(I)  Laser  simulator  panel 

(II)  Optical  contrast  simulator  panel 

(III)  Computer  remote  control  pane: 

(Iv)  Interconnection  patch  panel 

(v)  Missile  test  panel 

The  laser  and  optical  contrast  panels  each  contain  an  Input  selection 
switch  which  permits  Inpyt  voltages  to  originate  from  analog  computer 
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Figure  2.26  Dual  Glmbal  Hlrror  Calibration  - Up/Down 
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output,  manual  control  or  external  source.  When  manual  control  Is  selected 
three  potentiometers  on  each  panel  control  the  value  of  Input  voltages  to 
the  two  axes  for  each  mirror  and  for  zoom  and  spot  size  controller 
respectively  In  laser  and  optical  simulators.  In  manual  or  external  control 
voltage  Inputs  are  In  the  range  +10  volts. 

The  computer  remote  control  consists  of  a series  of  momentary  push- 
button switches  which  permit  the  transmission  of  logic  signals  to  the 
Cl-5000  logic  patchboard;  these  logic  signals  may  then  be  used  to  generate 
an  Interrupt  or  sense  line  Input  to  the  digital  computer.  Each  button  may 
be  lit  from  the  rear  under  computer  control  to  indicate  the  simulation  status. 

The  Interconnection  patch  panel  Is  a 37“3^  hole  Interchangeable 
patchboard  for  cross-connecting  various  elements.  The  holes  are  labelled 
horizontally  A to  Z,  AA  to  HM  and  vertically  1 to  3*»  (from  the  top  down). 
Included  In  the  control  console  are  sixteen  operational  amplifiers  fe*  use 
as  summers/ Inverters.  Inputs  and  outputs  for  these  appear  on  the  patchboard. 
Tables  IX  and  X Indicate  the  allocation  of  some  of  the  patchboard  connections. 

The  control  console  Includes  a missile  test  panel  which  cental*;  test 
point  outputs  and  control  Input  monitoring  points  for  a particular  set  of 
missile  hardware.  These  monitoring  and  Input  polntr.  are  connected  to  the 
patchboard  as  shown  In  Table  X. 
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TABLE  IX  PATCHBOARD  CONTENTS 


Connection  Element 

Patchboard 

Coordinates 

Connection  Element 

| 

Patchboard 

Coordinates 

Cl-5000  Trunklines 

Laser  Simulator 

Control  Panel 

T066 

A-l 

T067 

A-2 

Mirror  Pitch  (computer) 

C-ll* 

T068 

A-3 

Mtrrow  yaw  (computer) 

C-15  ! 

T069 

A-A 

Spot  size  (computer) 

C-16  ! 

T070 

A-5 

Mirror  pitch  (External) 

C-17  1 

T07I 

A-6 

Mirror  yaw  out 

C-18  : 

T072 

A-7 

Spot  size  out 

C-19  ! 

T073 

A-8  1 

Mirror  pitch  (monitor) 

C-ll 

T07^ 

A-9 

Mirror  yaw  (monitor) 

c-12  : 

T075 

A- 10 

Spot  size  (monitor) 

c-13 

TH76 

A-ll  :!  Ground 

C-20 

T077 

A-12 

T078 

A-13 

Optical  Contrast  Control 

T079 

A-l* 

Panel 

T080 

A-15 

T081 

A-16 

Mirror  pitch  (computer) 

C-2 1 

T082 

A-17 

Mirror  yaw  (computer) 

C-22 

T083 

A-18 

Zoom  (computer) 

C-23 

T084 

A-19 

Mirror  pitch  (external) 

G-14 

T085 

A-20 

Mirror  yaw  (external) 

G-15 

TO86 

A-21 

Zoom  (external) 

G-16 

T087 

A-22 

Mirror  pitch  out 

C-17 

TO88 

A-23 

Mirror  yaw  out 

C-18  1 

T089 

A-2  h 

Zoom  out 

C-19  1 

T090 

A-25 

Mirror  pitch  (monitor) 

c-2*  : 

T091 

A-26 

Mirror  yaw  (monitor) 

c-25 

T092 

A-27 

Zoom  (monitor) 

c-26 

T093 

A-28 

Ground 

G-20 

Cl -500  Trunklines 

Missile  Test  Panel 

T1<|0 

C-i 

Test  Point  12 

G-l 

Tll*1 

C-2 

Test  Point  13 

G-2 

Tl*2 

C-3 

Test  Point  1*» 

G-3 

T143 

C-4 

Test  Point  15 

Q-h 

TIM 

C-5 

Test  Point  16 

G-5 

J\kS 

C-6 

Test  Point  17 

G-6 

J\k6 

C-7 

Test  Point  26 

G-7 

TU7 

C-8 

Test  Point  27 

G-8 

TI50 

C”9 

Ground 

G— 10 

T151 

C-10 

Vane  1 

E-l 

TABLE  IX  PATCHBOARD  CONTENTS  cont 


Patchboard 

! i 

Patchboard 

Connection  Element 

Coordinates 

Connection  Element 

Coordinates 

Missile  Test  Panel 

i 

cont . 

i 

■ i 

i ] 

! 

Vane  2 
a 
a 

Vane  Ground 
Orlv  I 
Drfv  2 
Drlv  3 
Orlv 

Orlv  Ground 
Roll  Monitor 
Yaw  rate 
Pitch  rate 
Roll  gyro 
Seeker  yaw 
Seeker  pitch 
Launch  timer 
Forward  Loop. 
Forward  Loop" 
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SOFTWARE 


3.1  Mon Itor/Opera ting  System 

All  digital  and  hybrid  (except  special ly  wrl tten  'stand-alone')  pro- 
grams execute  on  the  Sigma  5 under  the  control  of  an  operating  system  pro- 
gram (also  called  a monitor  system).  The  monitor  currently  used  Is  an 
edcptatfon  of  a standard  Xerox  software  product  called  RBM  and  Is  known 
as  HRBM  (for  Hybrid  Real-Time  Batch  Monitor).  HRBM  handles  all  Input/ 
output  through  the  M10P,  handles  job  start  and  termination,  controls  the 
allocation  of  core  memory  and  use  of  disk  storage  space,  provides 
Interrupt  connect  lon/re l ease/processing  Including  loading  and  executing 
real-time  tasks.  A detailed  description  of  RBM  Is  given  In  References  6, 

7 and  8 and  of  HRBM  In  Reference  9.  An  overview  of  HRBM  will  be  given  here 
and  the  modifications  to  RBM  In  forming  HRBM  will  be  Indicated. 

Core  memory  Is  divided  Into  three  regions  described  as  follows  In 
order  of  Increasing  addresses: 

(I)  Core-resident  HRBM 

(it)  Background  area 

(til)  Foreground  area 

The  division  of  memory  not  used  by  (1)  Into  regions  (11)  and  (III)  Is 
determined  Initially  at  the  time  the  i.ionitor  system  Is  generated,  but  Is 
alterable  by  operator  key-ln  at  the  teletype  (RBM)  and  by  a program  control 
card  (HRBM  only);  If  the  foreground  area  is  to  be  truncated  but  contains 
a program  the  the  alteration  Is  not  accepted  by  the  monitor  and  Is  Ignored; 
a message  to  this  effect  Is  typed  on  the  operator  console  (teletype). 

Two  batch  Job  streams  are  permitted  In  HRBM  (one  In  RBM),  which 
execute  (not  simultaneously)  in  the  background  area.  The  two  streams  are 
labelled  'local'  and  'remote*  with  all  local  Jobs  having  priority  over 
remote  Jobs.  A card  reader  and  line  printer  Is  allocated  to  each  stream 
and  readers  and  printers  may  be  switched  by  operator  key-ln.  If  a remote 
Job  Is  executing  and  a local  Job  Is  entered  at  the  local  card  reader  then 
the  remote  Job  Is  suspended  (known  as  'checkpointed'  In  References  6,  7, 

8 and  9),  the  memory  contents  and  CPU  context  saved  on  a RAD  file,  and 

the  local  Jobs  processed.  The  remote  job  Is  restored  from  the  RAD  file  and 
execution  continued  after  all  local  Jobs  have  been  completed.  Since  there 
Is  no  protection  of  the  I/O  devices  which  can  be  accessed  by  both  local 
and  remote  jobs,,  care  must  be  exercised  to  avoid  a checkpointed  job  being 
Jeopardized  by  Interference  with  its  uncompleted  use  of  an  I/O  device.  I/O 
devices  In  this  category  are  the  magnetic  tape  d Ive  and  the  graphics 
display  unit. 

* 
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The  foreground  region  is  used  by  programs  which  respond  to  external 
Interrupts  and  thus  may  interrupt  execution  of  background  processing 
(which  may  be  local  or  remote).  Foreground  programs  may  occupy  any  part 
of  the  combined  background/foreground  areas;  when  a foreground  program 
uses  a part  of  the  background  region  an  execution  background  program  is 
checkpointed  as  described  above  but  on  a different  RAD  file.  Foreground 
programs  may  be  Initiated  from  the  background  job  streams,  or  from 
previously  compiled  and  loaded  RAD  files  by  operator  key-in.  All  com- 
pilers, assemblers  and  loaders  (or  link-editors)  execute  as  background 
Jobs. 


Card  reader  Input  and  line  printer  output  in  HRBM  is  handled  through 
what  are  called  'I/O  symbionts'.  This  means  that  the  card  inputs  and 
print  outputs  are  buffered  via  RAD  files;  all  cards  are  input  via  HRBM  to 
the  RAD  buffer  file  and  all  requests  to  read  a card  by  a background  pro- 
gram (foreground  programs  are  not  permitted  to  read  cards)  are  diverted  to 
read  the  next  available  card  image  from  the  RAD  buffer.  Similarly,  printer 
output  is  Intercepted  in  HRBM  and  stored  on  a RAD  file.  Output  to  a 
printer  thus  occurs  asynchronously  with  the  program  originating  the  output 
at  a rate  dependent  on  print  line  availability  and  printer  speed.  The 
symbiont  I/O  system  is  driven  by  internal  clock  3 which  interrupts  pro- 
cessing every  200  milliseconds  to  initiate  card  reading  and  line  printing 
on  both  readers  and  both  printers  if  required.  Once  Initiated,  reading 
and  printing  are  continued  by  t.  e t/Q  Interrupt  response  routines.  If, 
however,  the  clock  3 Interrupt  routine  finds  an  external  Interrupt  of 
priority  higher  than  10  active  then  no  Initiatin'1  action  is  taken  in  order 
to  minimize  interference  with  high  priority  foreground  jobs.  Local  back- 
ground and  foreground  printer  outputs  use  different  RAD  file  buffers  but 
share  the  same  printer. 

Most  components  of  RBM  use  re-entrant  coding  which  thus  precludes 
the  use  of  extra  register  blocks.  HRBM  uses  the  extra  three  register 
blocks  mentioned  in  Section  2.1.1  as  follows: 

(I)  Register  Block  1 is  assigned  tc  the  clock  3 Interrupt  routine 
where  re-entrancy  is  prevented  by  disabling  the  clock  3 
counter  equals  0 interrupt  until  Just  prior  to  exit  from  the 
Interrupt  routine. 

(II)  Register  Block  2 is  assigned  to  the  clock  k Interrupt  routine 
and  the  control  panel  Interrupt  . outine.  Clock  i*  Is  used  to 
maintain  time  of  day  and  the  date. 

(ill)  Register  Block  3 Is  assigned  to  the  DMA  channel  interrupt 
response  routine. 

Extensions  to  the  list  of  unsolicited  key- Ins  for  HRBH  are  Indicated 
tn  Table  XI.  Note  that  the  SS  key-ins  are  only  necessary  to  Intervene 
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TABLE  XI  HRBM  Operator  Key! ns 


Keyln  Characters* 

Effect 

MST  V 
B 

Permit  a foreground  (F)  or  background 
(8)  program  to  select  Master  mode  In 
the  CPU 

SLV  F 

Inhibit  a foreground  (F)  or  background 

_B_ 

(B)  program  from  selecting  Master  mode 

In  the  CPU 

SDS  V 

Exchange  the  designation  of  local  and 

R 

remote  printers  (P)  or  card  readers  (R) 

SS  >' 

initiate  foreground  (F)  or  background 

B 

(8)  output  symbiont  printing 

* Brackets  []  Indicate  a character  to  be  selected  which  must  be 
preceded  by  a blank. 


in  the  current  output  to  the  local  line  printer  in  order  to  switch  the  out- 
put from  foreground  to  background  or  vice  versa.  Under  normal  circumstances, 
ss  soon  as  the  RAO  output  buffer  being  printed  is  empty  the  symbiont  I/O 
control  routine  switches  to  the  other  output  buffer.  If  simultaneously  core- 
resident background  and  foreground  jobs  are  both  producing  printer  output 
this  may  result  in  the  two  sets  of  outputs  being  intermixed.  It  may  be 
necessary  to  direct  the  foreground  output  to  a RAD  file  and  then  print  the 
file  by  means  of  a separate  background  job. 

As  an  added  convenience  in  HRBM  the  SY,  FG  and  FMEM  keyins  of  RBM  have 
been  extended  to  become  control  commands  with  formats: 

I5Y 

!FG 

IFMEM  n 

and  are  Inserted  In  a job  deck  at  the  points  where  the  keyins  would  normally 
be  required.  Note  that  for  the  FMEM  command  the  foreground  area  should  not 
be  occupied  if  I *■  Is  being  allocated  to  background,  otherwise  the  command  is 
Ignored  with  an  operator  message  to  this  effect.  These  control  card 
commands  have  e:  actly  the  same  effect  as  the  corresponding  key  in  commands. 

A further  extension  included  in  HRBM  are  CAL  Instruction  responses  to 
meet  the  requirements  of  the  input/output  control  from/to  the  graphics 
display  terminal;  this  control  includes  central  connection  of  and  response 
to  external  Interrupts  6 and  7 for  character  receipt  and  transmission. 

Further  details  are  given  In  Section  3.5. 

3.2  Compilers  and  Assemblers 

The  system  software  contains  two  FORTRAN  compilers  and  an  assembler, 
briefly  described  as  follows: 

(I)  Extended  H level  FORTRAN  (referenced  by  file  name  FORTRANH  in 
the  SP  area  of  RAD  storage).  This  compiler  Is  described  In 
Reference,  10  and  Is  a superset  of  standard  FORTRAN  IV.  It  Is 
a one-pass  compiler  designed  to  be  compatible  with  IBM  360 
H level  FORTRAN;  however,  In-line  assembly  language  coding  may 
be  Introduced  if  required. 

(U)  Extended  FORTRAN  IV  (referenced  by  file  name  FORTRAN  In  the  SP 
area  of  RAD  storage)  Is  also  a superset  of  standard  FORTRAN  IV 
but  contains  many  more  extended  features  than  H FORTRAN. 

Add! t icnal 1y , it  is  designed  for  rapid  compilation  and  produces 
highly  efficient  compiled  code.  Reference  11  contains  details 
of  the  compiler  and  its  use.  This  FORTRAN  compiler  requires  a 
larger  amount  of  core  space  than  does  the  H level  compiler, 
and  also  makes  greater  demands  on  the  temporary  background  (BT) 
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RAD  files.  For  a program  which  contains  one  or  more  relatively 
large  subroutines  It  Is  necessary  to  increase  the  allocated 
size  of  the  XI  scratch  file  with  a card  wh'ch  precedes  the 
IFORTRAN  card  and  Is  punched  as  follows: 

1ALL0BT  (F!LE,X1) , (FSIZE,200) , (RSIZE,202) 

For  very  large  routines  It  Is  sometimes  necessary  to  insert  a 
similar  card  for  the  X2  file. 

(Ill)  Assembler  (referenced  by  file  name  MACRSYM  In  the  SP  area  of 
RAD  storage).  This  Is  a two  pass  assembler  v/Ith  extenstve 
macro  facilities.  Input  (i.e.  assembly  language  programs)  may 
be  In  the  form  of  source  records  (card  images)  or  of  compressed 
symbolic  records.  In  the  latter  case  an  updating  facility  Is 
provided  and  output  of  a compressed  symbolic  deek  may  be  speci- 
fied. Further  details  are  given  in  Reference  12.  The  language 
is  named  MACROSYMBOL;  a useful  feature  is  the  capability  of 
automatically  Including  within  a program  the  contents  of  a 
source  file  on  the  RAD  by  means  of  a SYSTEM  directive  and  a 
control  card  entry. 

3.3  Hybrid  Computer  Control  and  1/0  Routines 

A set  of  library  routines  is  provided  whtch  enables  a FORTRAN  or 
MACROSYMBOL  program  to  communicate  with  the  analog  computers  via  the  A/D 
Interfaces.  These  library  routines  are  formed  Into  a pre-complled  and 
linked  Public  Library  which  Is  stored  In  a RAD  file  named  HYBLIB  In  the  FP 
area.  To  use  this  Public  Library  a field  containing  (PUBLIB,  HYBLIB)  must 
be  Included  In  the  users  I0L0AD  card  (see  Reference  6,  Section  6).  The 
Public  Library  Is  capable  of  being  shared  by  foreground  and  background 
programs.  A description  of  the  method  of  use  of  these  routines  Is  given 
In  Reference  2. 

The  hybrid  library  routines  permit  communication  with  both  analog 
computers  (C I -500  and  Cl-5000)  and  provide  control  of  the  analog-digital 
Interfaces  and  the  high-speed  DMA  channel.  A brief  list  of  the  operations 
performed  Is  as  follows: 

I)  read  any  addressable  analog  component 

II)  set  potentiometers 

III)  read  analog  status  lines 

Iv)  read  the  digital  voltmeter 
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v)  set  the  analog  reference  dividers 
v I ) control  the  analog  modes 
vll)  send  analog  codes 
vl It)  write  MDAC's 
lx)  read  ADC's 

x)  control  the  real-time  clocks 

xi ) read  sense  lines  (discretes) 
xll)  write  control  lines  (discretes) 
mill)  control  external  Interrupts 

For  real-time  operation  the  concept  of  "patterned"  execution  Is  Imple- 
mented. In  this  usage,  patterns  of  hybrid  I/O  and  control  operations  are 
defined  prior  to  real-time  execution  as  a series  of  computer  Instructions  In 
the  user  program  temp  stack.  Then,  during  real-time  operation,  the 
pattern  Is  invoked  which  causes  execution  of  the  pre-stored  stack  of 
Instructions  resulting  In  minimum  overhead  time  during  real-time  operation. 

The  hybrid  computer  library  routines  use  the  CAL3  and  CAL**  Instructions 
and  therefore  an  assembly  language  user  program  should  not  use  these 
Instructions  for  an y other  purpose.  Hybrid  library  routines  contain  exten- 
sive error  checking  and  diagnostic  outputs  to  crver  Illegal  use  of  their 
entry  points.  Construction  of  the  library  permits  the  concurrent  use  of 
the  routines  by  various  foreground  and  background  programs.  However, 
connection  of  interrupts  already  assigned  to  a program  Is  prevented,  as  Is 
the  arming  and  enabling  of  Interrupts  which  have  not  been  connected  through 
the  Interrupt  connecting  routine. 

The  task  control  block  (TCB)  for  each  connected  interrupt  Is  also 
allocated  to  the  user  program  temp  stack.  Thus,  the  use  of  patterns  and 
Interrupts  should  be  considered  when  allocating  a user  program  temp  stack 
Size. 


3.4  CSSL 

CSSl  (Continuous  Systems  Simulation  Language)  Is  a digital  computer 
problem-oriented  language  for  the  representation  of  continuous  dynamic 
systems  which  can  be  modeled  by  sets  of  ordinary  differential  equations, 
language  definitions  and  methods  of  use  are  given  In  Reference  13.  A 
program  formed  from  CSSl  statements  Is  converted  by  the  CSSL  translator 
Into  a FORTRAN  program  which  Is  then  compiled,  loaded  and  executed  In  the 
usual  manner. 


CSSL  permits  both  equation  based  and  block  oriented  representation 
statements  and  is  thus  suitable  for  a wide  range  of  users. 

The  language  translator  implemented  for  the  Sigma  5 requires  three 
scratch  files  in  the  background  temporary  area  of  the  RAr;.  These  files 
are  X4,  X5  and  X6  with  blocked  format,  record  size  128  •^•-ds  and  file  size 
400  records  for  an  average  program.  A large  program  v require  Increased 
file  sizes.  The  translated  output  is  a series  of  FORUvu*  statements  and 
the  output  device  is  usually  a magnetic  tape;  the  translated  program  is 
also  a foreground  program  with  the  first  subroutine  connected  to  interrupt 
14  (Interrupt  address  6Ejg).  However,  if  real-time  execution  is  not  re- 
quired the  first  three  statements  of  the  output  may  be  deleted,  in  which 
case  the  FORTRAN  output  may  be  executed  as  a background  job,  as  determined 
by  the  10L0AD  card. 

The  contrel  card  set-up  for  a CSSL  program  is  as  follows: 

! JOB  082274, ACJ 
!REW  9TA80 

’ALLOBT  (FI LE,X4), (FORMAT, B),(RSIZE, 128), (FSIZE, 400) 

1ALL0BT  (FILE, X5), (FORMAT, B),(RSIZE, 128), (FSIZE, 400) 

JALL0BT  (FILE.X6) , (FORMAT, B) , (RSIZE, 128) , (FSIZE, 400) 

1CSSL 


CSSL  Source  Cards 


• 

1WE0F  9TA80 
1REW  9TA80 
IPREC  9TA80.3 
! ASSIGN  (H:SI,9TA80) 

1 FORTRAN  SI ,G0,BC,NS,LS 
! 0L0A0  GO , (TEMP , 400) , (LI B , USER , SYSTEM) 

IROV 

1FIN 

Note  that  the  IPREC  card  skips  the  first  3 statements  to  avoid 
connecting  the  program  to  an  Interrupt.  If  foreground  operation  is  required 
the  (FORE)  field  Is  required  on  the  I0L0A0  and  an  1FG  card  is  required 
before  the  1R0V  card. 

3*5  Graphics  Software  i 

As  mentioned  In  Table  I,  the  graphics  terminal  Is  a Tektronix  Model 
4002A.  The  terminal  comprises  a storage  tube  display,  keyboard  input, 
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hardcopter  and  a "Joystick"  driven  cross-hair  cursor.  The  40Q2A  Is  con- 
nected to  the  Sigma  5 by  means  of  a Xerox  7611  character  oriented 
communications  controller  which  transmits  and  receives  ASCII  8 bit 
characters.  External  Interrupts  6 and  7 (interrupt  locations  66^  and 
67jz)  are  used  by  the  character  receive  and  transmit  circuitry  respectively 
to  Indicate  completion  status  to  the  Sigma  5.  Character  transmission  rate 
is  currently  set  at  9600  Baud;  the  7611  has  a capacity  of  64  lines  of 
which  only  one  (line  no.  0)  is  used  by  the  4002A. 

Standard  software  which  facilitates  the  use  of  the  graphics  terminal 
Is  provided  in  three  packages,  which  are  briefly  described  as  follows: 

I)  Hardware  I/O  Control  and  Utility  Routines. 

This  package  comprises  a set  of  FORTRAN  callable  assembly 
language  routines  (MACROSYMBOL)  which  perform  the  software 
Interfacing  function  between  the  Sigma  5 and  the  4002A.  Entry 
point  names,  parameter  definitions  and  operating  functions  may 
be  described  as  follows: 

CALL  ST7611  - this  entry  point  Initializes  the  software  package 
by  resetting  the  7611  and  checking  that  it  is  switched  on. 
Interrupts  6 and  7 are  connected  centrally  and  various  counters 
Initialized.  This  entry  point  should  be  called  once  before  use 
of  the  following  entry  points. 

CALL  WT761 1(1)  - transmits  one  ASCII  character  contained  in  the 
rightmost  byte  (l.e.  byte  3)  of  I to  the  display  screen.  The 
character  may  be  of  7 bit  form  without  the  parity  bit;  this  is 
calculated  and  inserted  by  the  routine.  Note  that  on  return 
from  this  call  the  character  has  not  necessarily  been  trans- 
mitted to  the  screen  since  the  output  Is  buffered. 

CALL  RO761 1 (N, I CHAR, IWT)  - reads  N characters  from  the  terminal 
and  stores  them  In  consecutive  bytes  In  array  ICHAR.  Parameter 
IWT  Is  used  to  Indicate  whether  it  is  desired  to  return  control 
to  the  calling  routine  before  or  after  the  read  operation  Is 
complete.  If  IWT  Is  non-zero  control  is  not  returned  until  the 
requested  number  of  characters  have  been  received,  and  If  IWT 
Is  zero  control  Is  returned  Immediately  after  the  Input  has 
been  Initiated.  In  the  latter  case  the  NSTCK  entry  point  may 
be  used  to  determine  the  status  of  the  I/O  operation. 

CALL  HT76I!  - Is  the  terminating  call  to  the  software  package 
which  stops  I/O  activity  and  disconnects  the  external  Interrupts. 
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The  remaining  entry  points  are  used  as  FORTRAN  function  calls; 

NEBCD(1,J)  - returns  the  EBCDIC  character  corresponding  to  the 
8 bit  ASCII  character  contained  In  the  J'th  byte  of  array  I 
counting  from  1 at  the  left  of  the  first  word.  Thus,  characters 
read  by  the  RD76II  entry  point  into  I may  be  conveniently  con- 
verted Into  the  SIGMA  5 internal  code. 

NASCII(I,J)  - returns  the  ASCII  character  corresponding  to  the 
8 bit  EBCDIC  character  contained  in  the  J'th  byte  of  array  I 
counting  from  1 at  the  left  of  the  first  word. 

NDIG(I,J,K)  - when  K is  positive  the  j'th  byte  taken  from  array 
1 Is  converted  to  the  equivalent  binary  integer  If  it  is  an 
ASCII  numeric  character  but  otherwise  returns  a zero.  If  K is 
zero  the  least  significant  5 bits  of  the  J'th  character  are 
returned  as  a binary  integer  and  if  K is  negative  the  least 
significant  7 bits  are  returned. 

NSTCK(l)  - is  used  to  check  the  I/O  status  where  I is  a dummy 
parameter.  Values  returned  are  integers  with  the  following 
meanings: 

0 The  software  is  uninitialized 

1 Software  has  been  Initialized  but  no  I/O  Is  in  progress 

3 Characters  are  being  output  to  the  terminal 

? 5 Input  from  the  terminal  has  been  requested  but  the  required 

number  of  characters  have  not  been  received 

7 Both  input  from  and  output  to  the  terminal  are  currently 
In  progress 

Any  other  code  means  the  same  as  zero. 

II)  Tektronix  Terminal  Control  System  - 4002A 

This  package  consists  of  a set  of  FORTRAN  routines  which  faci- 
litate the  generating  of  displays  and  provides  convenient 
control  of  the  screen  contents.  Reference  14  contains  full 
details  of  the  operation  of  use  of  this  package. 

The  package  contains  routines  for  direct  and  relative  vector 
drawing,  controlling  and  positioning  alphanumerics  including 


tabs  and  margins  and  provides  support  for  the  scratchpad  and 
edit  modes.  The  concept  of  virtual  graphics  Is  also  Implemented 
which  permits  the  generation  of  displays  on  a virtual  screen 
larger  than  the  actual  screen.  Portions  of  the  virtual  display 
may  then  be  viewed  through  one  or  more  screen  "windows"  which 
may  rotate  or  scale  the  actual  display  relative  to  the  virtual 
window. 

ill)  Tektronix  Advanced  Graphing  II 

This  package  consists  of  a set  of  FORTRAN  routines  which  are 
useful  In  generating  graphs.  A.  larger  number  of  options  are 
included  In  an  attempt  to  provide  a flexible  means  of  drawing 
all  the  types  of  graph  normally  encountered.  Full  details  are 
provided  In  Reference  15. 

All  the  software  described  in  this  subsection  Is  stored  on  RAD  files 
In  the  D6  area  In  source  and  relocatable  binary  form.  File  names  are  given 
In  the  following  table: 


Software 

File 

Names 

Package 

1 

1 Source 

t 

Relocatable  Binary 

COC  I/O  Control 

1 MOCPACK 

MOCA&2BO 

Terminal  Control  System 

TCSGRAP 

GRAPHICS 

Advanced  Graphing  II 

AG2FILE 

AG2FiLB0 

TABLE  XII  GRAPHICS  SOFTWARE  FI IF  STORAGE 


Software  under  this  heading  consists  of  a background  processor  which 
performs  various  editing  and  transmission  functions  on  any  symbolic  source 
card  Images.  The  processor  Is  written  In  MACROSYMBOL  assembly  language 
and  Is  held  In  the  system  programs  area  of  the  RAD  In  core  Image  format 
In  a file  named  UPDATE.  It  may  thus  be  loaded  Into  core  memory  for  exe- 
cution by  an  IUPDATE  control  card. 

The  operations  of  this  editing  processor  are  controlled  by  a series 
of  commands  provided  by  the  user  and  read  by  UPDATE  from  the  C device. 
Input  and  output  source  records  may  be  In  compressed  or  uncompressed  form. 
The  compressed  symbolic  form  Is  that  used  by  the  HACROSYHBOL  assembler. 
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Source  files  are  considered  conceptually  to  consist  of  a set  of  'decks', 
the  delimiters  of  a deck  being  records  containing  END  (except  for  the  first 
deck  in  a file,  where  the  first  record  is  the  initial  delimiter).  Thus, 
UPDATE  is  designed  to  handle  editing  of  source  programs  which  may  contain 
FORTRAN  and  MACROSYMBOL  decks  Intermixed.  However,  because  of  the  different 
requirements  for  END  cards  in  FORTRAN  and  MACROSYMBOL  a command  is  provided 
to  indicate  which  type  of  END  card  is  used  as  a delimiter  for  the  current 
deck. 

UPDATE  control  commands  per?  «t  the  following  operations: 

. copy  decks  from  an  input  ' ■ an  output  device 

. during  a deck  copying  operation  insert  and  delete  records  in  the 
deck 

. control  the  selection  of  input,  output,  command  or  listing  devices 
during  deck  transferral  and  thus  permit  the  deletion  or  addition 
of  decks  to  a file 

. convert  source  code  from  IBM  026  code  to  the  029  code  if  required 

. write  end  of  file  on  the  output  device 

. rewind  selected  devices 

* 

. assign  DCB's  corresponding  to  operational  labels  to  RAD  files  or 
hardware  devices 

The  control  commands  to  perform  these  operations  are  described  as 
follows.  Note  that  each  control  command  format  Is  :X  followed  by  a single 
blank  where  X Is  one  of  the  coomands  described  below.  Punching  always 
starts  in  card  column  1. 

I)  :C  xx 

Specifies  the  cotrmand  input  device  where  xx  Is  its  operational 
label  and  may  be  C,  0C,  Cl  or  SI.  The  default  Is  C 

II)  :l  xx 

Specifies  the  ’nput  device  where  xx  Is  Its  operational  label 
and  may  be  Si  or  Cl.  The  default  Is  St  and  If  Cl  Is  specified 
the  source  must  be  in  standard  MACR0SYM60L  compressed  symbolic 
form. 
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m)  :o  xx 

Specifies  the  output  device  where  xx  Is  Its  operational  label 
and  may  be  SO,  CO  or  DO.  DO  Is  delected  output.  In  which  case 
the  decks  being  copied  are  deleted  until  the  output  device  Is 
changed.  CO  produces  standard  compressed  symbolic  output. 

The  default  Is  SO. 

Iv)  :L  xx 

Specifies  whether  or  not  the  output  file  listing  Is  required 
where  xx  Is  LO  for  listing  and  NO  for  no  listing.  The  default 
Is  NO.  The  listing  will  appear  on  the  system  LO  device  and 
will  have  the  Input  source  line  numbers  with  all  Insertions 
flagged. 

v)  :A  xx,yy,nnnnnnnn 

Assigns  the  DCB  corresponding  to  operational  label  xx  tc  device 
yy  (TV  for  TYA01,  CR  for  CRA03,  9T  for  9TA80)  or  to  RAD  area  yy 
(SP,  FP,  BP,  BT,  Dl,  ...,DF).  If  the  assignment  is  to  a RAD 
file  area  nnnnnnnn  specifies  the  file  name  but  if  yy  Is  a RAD 
area  and  the  file  name  Is  omitted  the  DCB  is  assigned  to  the 
entire  RAO  area,  xx  may  be  SI , Cl , OC,  C,  SO  or  CO. 

vl)  :FOR 

Specifies  that  decks  are  to  be  delimited  by  FORTRAN  END  cards. 
The  default  Is  MACROSYNBOL. 

vll)  :MAC 

Specifies  that  decks  are  to  be  delimited  by  NACR0SYH80L  END 
cards. 

vl 1 1 ) :E0F 

Writes  an  end  of  file  mark  on  the  output  device, 
lx)  :REV  xx 

Rewinds  the  device  assigned  to  operational  label  xx  where  xx 
may  be  St,  Cl,  SO  or  CO.  If  xx  is  unspecified  the  output 
device  will  be  rewound. 

x)  :B 

Specifies  that  026  to  029  conversion  "switch"  Is  to  be  toggled. 
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The  default  Is  no  conversion.  The  first  time  this  card  Is 
encountered  conversion  of  the  Input  begins  and  the  next  time 
a card  of  this  type  is  encountered  the  conversion  Is  suppressed, 
and  so  on. 

xi)  :END 

Returns  control  from  UPDATE  to  the  monitor. 


The  remaining  control  cards  refer  to  the  control  of  copying  decks 
within  a file.  These  cards  are  indicated  by  a + or  - in  column  1 of  each 
card. 

a)  +END  m 


Specifies  that  decks  are  copied 
vices  until  m decks  (i.e.  m END 
missing  m defaults  to  the  value 

b)  +E0F  m 


from  the  input  to  the  output  de- 
cards) have  been  copied.  A 
1. 


Specifies  that  decks  are  copied  from  the  input  to  the  output  de- 
vices until  m end  of  file  marks  have  been  copied.  A missing  m 
defaults  to  the  value  1;  the  file  marks  are  copied  to  the  output 
file. 


c)  -EOF  m 

Identical  to  b)  above  except  that  the  end  of  file  marks  are  not 
copied  to  the  output  file. 

d)  -Hn.n 

Specifies  that  line  numbers  m through  n In  the  current  deck  are 
to  be  deleted  during  copying  to  the  output  file. 

e)  4® 

Specifies  that  the  current  deck  Is  to  be  copied  to  the  output 
file  up  to  and  Including  line  number  ra. 

Any  record  read  from  the  command  device  which  cannot  be  Inter- 
preted as  an  UPDATE  command  is  written  to  the  output  device  at 
the  current  line  position.  For  this  reason  Insertion  and  Jeletlon 
commands  d)  and  e)  must  appear  In  the  set  of  controt  commands 
grouped  for  each  deck  In  ascending  numerical  line  number  order. 
Insertion  records  should  appear  immediately  following  the  type 
d)  or  e)  record  to  which  they  refer. 
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EXAMPLE  HYBRID  SIMULATION  PROBLEM 


As  a typical  example  of  the  application  of  the  analog/digital  hybrid 
computer  described  In  the  foregoing  sections  of  this  report,  a simulation 
of.  a simplified  missile  Intercepting  a target  In  two  dtmenslonal  space  Is 
given  In  this  section.  Motion  Is  assumed  to  occur  In  a vertical  plane, 
thus  requiring  two  translational  and  one  rotational  degrees  of  freedom  to 
describe  the  motion  of  the  missile.  The  target  is  assumed  to  be  stationary. 

The  simplified  mathematical  model  of  the  missile  syste'  ' described, 
followed  by  a discussion  of  the  partitioning  of  the  problem  ».'<  analog  and 
digital  sections.  The  design  of  the  digital  program  and  Its  I r» *■ ' action 
with  the  analog  computer  Is  described  In  detail  and  finally  a set  of 
results  Is  Included. 

4. 1 Mathematical  Model  Description 

4. 1.1  Vehicle  Geometry,  Mass  and  Inertia 

Reference  length  .584  ft 

Reference  area  .267  ft2 

Initial  mass  Mq  4.5942  slugs 

2 

Initial  pitch  axis  moment  of  Inertia  lyyo  6.9  slugs  ft 

* 

Mass  variation 

dm  * -Thrust  0 < t < t 

dt  g SI  burn  (4.1) 

- 0 t > t 

burn 

where  g Is  gravity  acceleration  and  SI  Is  thu  engine  specific 

impulse. 

Inertia  variation 


dm 

It 


(4.2) 


4.1.2  Engine  Thrust 

Engine  thrust  Is  3000  lb  along  the  missile  centerline  for  a 
duration  of  3 seconds  (l.e.  ^ ■ 3 secs). 

Engine  specific  Impulse  SI  ■ 250  secs. 


4.1.3  Axis  Systems  and  Reference  Frames 

Reference  frames  used  In  the  simulation  include  an  inertial  frame 
XOZ  with  the  origin  0 at  the  point  on  th  sarth  surface  given  by  the 
launch  point;  a frame  fixed  in  the  body  xoz  with  origin  at  the  missile 
CG  and  a frame  fixed  in  the  space  referenced  seeker.  The  inertial  X axis 
Is  horizontal  pointing  towards  the  target  and  the  inertial  Z axis  is 
vertically  downwards. 

The  angle  between  the  body  axes  and  inertial  axes  and  inertial  axes 
Is  measured  positively  when  Z is  rotated  towards  X.  The  angle  between 
the  seeker  frame  x axis  and  the  body  x axis- is  6^  are  both  zero  all  che 
reference  frames  coinc’de. 

4.1.4  Aerodynamic  Forces  and  Moments 


Aerodynamic  forces  and  moments  are  simplified  by  assuming  linearity 
with  angle  of  attack  x,  and  control  vane  angle  6,  and  are  expressed  In  the 
form  of  tabulated  functions.  The  functions  required  are  axtal  force 
coefficient  C^,  normal  force  derivative  C^,  pitching  moment  derivative 
Cmx>  pitch  damping  derivative  CmQ,  control  vane  force  derivative  0^,$. 

All  these  are  functions  of  the  Mach  number,  and  given  in  Table  XIII.  In 
this  table  all  derivatives  are  per  radian,  except  for  CmQ  which  is  per 
radtan/sec.  Moment  derivatives  are  referred  to  the  missrle  CG. 

The  variation  of  air  mass  density  and  speed  of  sound  with  altitude 
Is  approximated  by  the  following  expression: 

Table  XIII  Aerodynamic  Functions 


cx 

CN 

% 

.305 

8.94 

-7.16 

-348 

1.6125 

.305 

14.9 

-10.3 

-405 

2.15 

.54 

16.6 

- 9.78 

-410 

2.15 

.755 

14.75 

- 6.94 

-406 

2.29 

.832 

12.89 

- 2.29 

-336 

2.25 

where  f Is  the  air  density,  V Is  the  local  speed  of  sound  and  Z Is  the 
missile  Z coordinate  which  Is~norrnal ly  negative  (for  positive  altitude). 

Aerodynamic  forces  In  the  body  frame  X and  Z directions  are: 


FXA  e "9s  c (4.5) 

FZA  - -qS(CN/+  CN66)  (4.6) 

2 

where  q Is  the  aerodynamic  pressure  iPV  and  6 Is  the  control  vane  angle, 
positive  trailing  edge  down.  The  pitching  moment  about  the  CG  Is  given 
by 


My 


ciSctC^- 


C_  .6 


(4.7) 


where  Q.  is  the  missile  pitch  rate. 

4.1 .5  Translation  and  Rotation  Equations 

The  total  force  components  acting  on  the  missile  are: 


FXB  “ fXA  ' m9sIn6E  + Ftx  (4.8) 

FZB  " FZA  + m9cos6E  + Ftz  . (4.9) 

where  F~x  and  F_z  are  thrust  components  along  the  missile  x and  z axes 
respectively.  Relative  to  the  inert’al  frame,  the  translational  equations 
of  motion  are: 


x e PXB  cos6E  + fZB  s!n°E 
m 

Z * F^g  sineE  + FZg  cosO^ 
m 


m 


M 


o 


dm 

dt 


and  the  magnitude  of  the  missile  velocity  Is 


(4.10) 

(4.11) 


V - /x2  + z2 

Rotation:'1  equations  of  motion  are 
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(4.12) 


(4.13) 


‘YYo  ^ t 


(4.14) 

(4.15) 


The  aerodynamic  angle  of  attacks  may  be  defined  as 
W 

tan  - WB 


(4.16) 


where  Ug  and  Wg  are  velocity  components  along  the  x and  z body  axes. 
Alternatively,  the  angle  of  attack  may  be  Introduced  as  a state  equation 
given  by 


cos*  - Fyg  sin* 


(4.17) 


Equation  (4.17)  has  some  advantages  in  a hybrid  application  over  (4.16)  and 
will  be  used  In  the  simulation  program. 

4.1.6  Guidance  and  Control 

The  target  Is  assumed  to  be  designated  by  a pulsed  laser  operating 
at  20  pulses  per  second.  A laser  detector  is  mounted  on  a glmballed 
platform  (one  degree  of  freedom)  in  the  nose  of  the  missile;  the  detector 
receives  the  laser  pulses  reflected  by  the  target,  and  generates  guidance 
signals  from  them  using  a form  of  propo;tional  navigation. 

The  pulsed  laser  effect  Is  equivalent  to  providing  a sampled  data 
seeker  as  Illustrated  in  block  diagram  form  in  Figure  4,1. 


Line  of  sight)  j* 

Target  . 

to  target  In 

y 

Position 

j gyro  frame 

s. 

Limit 

12° 


Figure  4.1  Seeker  Block  Diagram 
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The  line  of  sight  from  missile  to  target  Is  defined  by  a line  of 
sight  angle  6...  In  the  reference  frame  measured  positively  In  the  same 
direction  as  These  angles  are  Illustrated  in  Figure  4. 2.  If  the 

missile  position  Is  {Xm,  Zm)  and  the  target  position  Is  (Xy»  Zy)  then 


’LOS 


tan 


-’Hzt-zm) 


and,  In  the  seeker  frame, 


3L0$ 


(4.18) 


(4.19) 


X Is  converted  to  an  error  signal  by  the  llmlttng  process  Illustrated  In 
Figure  4.1  and  Is  then  multiplied  by  Ks,  the  effective  gyro  precession 
rate,  to  give  a»slgnal  X* which  drives  the  gyro  platform  and  the  rear- 
mounted  missile  control  fins  to  reduce  X to  zero  and  to  align  the  missile 
along  the  seeker  axis.  The  gyro  equation  Is  simplified  to 


(4.20) 


and  the  autopilot  block  diagram  In  Laplace  plane  form  is  given  In  Figure 
4.2.  Constant  terms  In  Figures  4.1  and  4.2  have  the  following  values: 


KA  - 2.  Kg  - .1,  Kfl  - 60,  K„  - .3703,  1^  - 3 
t,  » .025,  x2  “ -005,  t j « .01,  * 0, 

L ■ .6,  a>n  ■ 60 

The  0 Input  to  the  autopilot  In  Figure  4.3  I?  a pitch  damning  loop 
where'n  8 Ts  differentiated  and  smoothed  to  approximate  the  body  pitch 
rate  U and  the  lead-lag  network  filters  the  signal  to  avoid  exciting  the 
missile  natural  frequency  in  pitch. 

4.1.7  Target  Intercept 

Intercept  occurs  during  an  Integration  step  In  which  either  of  the 
missile  position  coordinates  exceeds  the  corresponding  target  coordinate. 
Thus,  at  the  end  of  each  integration  step  the  missile  position  is  extra- 
polated over  the  next  step,  using  the  current  velocity  to  test  whether 
Intercept  will  occur,  l.e.  If 

X ♦ X At  > XT  (4.21) 

m m i 

or  Zm  + Zra  At  > ZT  (4.22) 


% 


4-5 


If  either  of  the  above  Inequalities  Is  satisfied  the  simulation  Is  termin- 
ated at  the  current  step  (time  T) . The  miss  distance  over  the  next  time 
Interval  Is 


Rffllss 


(T-t)}2  * (ZT  - Zm 


(4.23) 


where  t Is  the  time  within  the  next  Interval  T+6t.  f^z  Is  a minimum  at 
the  Point  of  closest  approach  of  missile  to  target;  differentiating  (4.23) 
and  equating  to  zero  gives  a value  of  t at  intercept  of 


<XT  - 3.  * VI  *_  2m  <*T  - 2„  * 2J>  T 


where  the  subscript  T Implies  evaluation  at  time  T.  Substituting  t given 
by  (4.24)  Into  (4.23)  gives  the  approximate  miss  distance,  neglecting 
acceleration  over  the  last  partial  time  step. 

4.2  Partitioning  of  the  Problem 

Partitioning  of  the  mathematical  model  between  the  digital  and  analog 
computers  Is  constrained  by  the  requirement  that  the  simulation  execute 
In  real-time.  This  constraint  Implies  that  If  all  the  equations  are 
Implemented  on  the  digital  computer  all  the  12  states  (4  translational,  2 
rotational  and  6 In  the  seeker,  autopilot  and  actuator  combination)  must  be 
Integrated  In  real-time.  On  the  Sigma  5 described  in  the  foregoing  sections 
this  would  only  be  possible  using  large  Integration  steps  and  a simple 
Integration  process.  However,  the  time  constants  In  the  autopilot  are  such 
as  to  require  an  Integration  step  of  the  order  of  a millisecond.  Thus, 

It  Is  necessary  to  allocate  part  of  the  mathematical  model  for  simulation 
on  the  analog  computer  to  take  advantage  of  Its  high  speed  parallel  compu- 
tation properties. 

The  parts  of  the  mathematical  model  most  appropriate  for  assignment 
to  the  analog  computer  are  those  containing  the  highest  natural  frequencies 
and  short  time  constants.  The  actuator,  with  wn  **  60  rad/sec  Is  an  obvious 
candidate,  as  are  the  functional  blocks  In  Figure  4.2  which  operate  on  the 
gyro  angle  0^  and  line  of  sight  rate  Inputs  to  the  autopilot.  The  time 
constants  of  25,  5 and  10  milliseconds  there  represent  break  frequencies  of 
40,  200  and  100  rad/sec  respectively. 

The  seeker  sample  and  hold  process,  however,  Is  suitable  for  Imple- 
mentation on  the  digital  computer  since  this  has  a period  of  50  milliseconds 
(20  pulses/second)  and  a 0AC  output  Is  Ideally  suited  to  providing  the 


sampled  values.  An  objection  to  this  can  arise  In  the  sampling  of  the  gyro 
angle  6 which  will  be  sma'l  over  most  of  the  flight  but  must  be  scaled  to 
the  maximum  limits  on  0$  (normally  around  20°).  To  overcome  this  dis- 
advantage the  gyro  Is  Implemented  on  both  digital  and  analog  computers.  On 
the  analog  computer  0 Is  calculated  according  to  equation  (4.21)  and 
used  as  an  Input  to  tne  autopilot.  On  the  digital  computer  the  gyro  angle 
Is  calculated  relative  to  inertial  axes  to  give 

6S|  * At  Z X*  (4.25) 

where  At  is  the  50  millisecond  sample-hold  Interval  and  X Is  the  line  of 
sight  sample  converted  to  a rate.  From  equation  (4.25)  the  line  of  sight 
angle  In  the  seeker  frame  becomes 


X 


LOS 


- e 


si 


(4.26) 


The  missile  translation  equations  are  assigned  to  the  digital  com- 
puter because  of  the  low  frequencies  contained  In  the  velocity  and  dis- 
placements and  because  of  the  need  for  a high  dynamic  range  required  to 
give  an  accurate  miss  distance  at  the  end  of  a flight.  There  remains  to 
be  assigned  then  the  missile  rotational  equations  for  0 and  0£.  Q wl 1 1 
respond  approximately  as  a damped  second  order  system  with  a natural  fre- 
quency given  by  (see  reference  lb): 


(4.27) 

(4.28) 


Assuming  conditions  of  V » 1000  ft/sec  at  an  altitude  of  1000  ft  and  sub- 
stituting the  appropriate  values  after  engine  burnout  Into  equations  (4.27) 
and  (4.28)  gives  u.  of  17.1  rad/sec  which  Is  Just  less  than  3 Hz.  A 
simple  Integrat Ion '‘process  using  one  derivative  evaluation  per  step  would 
require  an  Integration  Interval  of  less  than  10  milliseconds  which  would 
probably  permtt  Q and  0p  to  be  assigned  to  the  digital  computer  In  this 
simple  3 degree  of  freedom  problem.  For  a full  dof  case,  however,  the 
time  requirement  would  prevent  this,  so  to  make  the  present  problem  a 
realistic  example  the  equations  for  Q,  0^  and  0S  are  assigned  to  the  analog 
computer.  In  addition,  </.  Is  evaluated  according  to  equation  (4.17)  on  the 
analog  computer  and./.,  A and  Q are  used  as  analog  DAC  multipliers  in  the 
evaluation  of  the  aerodynamic  moment  Input  to  equation  (4.13)  with  the 
aerodynamic  derivatives  Cm  , Cm  and  C.,  obtained  by  table  look-up  on  the 
digital  computer.  Evaluating  Hy  In  th I ft  manner  reduces  the  "staircase" 
effect  of  the  0AC  outputs. 


I 


Four  translation  states  are  Integrated  on  the  digital  computer  and 
all  body  force  and  moment  terms  are  evaluated  there.  The  digital  Inte- 
gration method  Is  a 2 point  Adams-Bashforth  process  with  projection  of 
derivatives  from  the  previous  time  frame,  which  may  be  represented  by  the 
following  equations  for  the  X Inertial  direction: 


*1+1 


:c,  +7<*l  + x'i-l! 

(4.2S) 

Cm  * y <«l  * *M> 

(4.30) 

+ 5-  (xr  + x ) 

1 2 ci+l 

' length  Is  h,  1 refers  to  the  current 

(4.31) 

time  and  the  suffix 

M+l 


C Implies  a "corrected"  value. 


Termination  of  the  simulation  and  the  calculation  of  miss  distance 
employs  equations  (4.21)  through  (4.23). 

A block  diagram  of  the  partitioned  problem  Is  given  In  Figure  4.3. 

4.3  Design  of  the  Digital  Program 

The  design  of  the  digital  section  of  the  hybrid  program  centers  around 
the  timing  and  Interrupt  sequencing.  The  real-time  section  of  the  simu- 
lation Is  driven  by  a series  of  clock  Interrupts  which  synchronize  the 
program  with  real-time  and  provide  control  of  the  digital  Integration 
process.  A convenient  reference  cycle  within  the  simulation  Is  the  50 
millisecond  seeker  sample-hold  Interval  and  the  digital  Integration 
Interval  Is  chosen  to  be  an  Integral  multiple  of  this  overall  time  frame. 

A digital  Integration  step  of  25  milliseconds  Is  adequate  for  the 
translation  equations  which  Is  also  the  Interval  for  output  of  the  aero- 
dynamic moment  update  terms  and  for  the  output  to  the  angle  of  attack 
equation.  A diagram  Illustrating  time  usage  during  the  real-time  section 
of  the  program  Is  given  In  Figure  4.4.  All  the  program  sections  are 
written  in  FORTRAN. 


In  addition  to  the  real-time  section  the  program  contains  Initiali- 
zation, pre-real  time  and  post-real  time  sections.  The  initialization 
section  Is  the  Main  Program  and  Is  executed  only  once  at  the  time  the 
program  Is  loaded  Into  the  foreground  region  of  memory.  Initialization 
consists  of  connecting  the  Sigma  5 external  Interrupts  to  various  routines 
In  the  program  and  defining  the  hybrid  1/0  patterns.  Thereafter,  for  all 
subsequent  simulation  runs  the  program  Is  driven  by  the  external  Interrupts. 


* 


( 
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rigure  4.3  Partitioned  Problem 
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Integration 
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Routine  SCHEDULE  (plus  time  unused) 
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Seeker  sample-hold  calculation  and  DAC  output 
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Figure  4.4  Frame  Time  Usage 
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The  pre-real  time  section  consists  of  subroutine  STARTRUN  whtch  Is 
connected  to  Interrupt  14.  This  subroutine  Initializes  various  digital 
and  analog  variables  and  reads  various  run  control  parameters  from  a 
RAO  file.  Run  control  parameters  are  read  from  file  SAMPLE  In  the  07 
RAD  area  and  are  punched  as  for  a NAMELIST  Input  (see  references  10  and 
13)  as  given  In  Table  XIV.  Variables  which  are  Included  In  Table  X V, 
but  which  are  not  read  from  file  SAMPLE  have  the  default  values  Indicated 
In  the  table.  At  the  <*nd  of  subroutine  STARTRUN  an  Initial  Integration 
step  Is  taken  by  triggering  the  appropriate  interrupt,  the  real-time 
clock  Is  then  set  running  and  the  analog  computer  is  set  to  the  compute 
mode  which  thus  commences  the  simulation. 

Subroutine  CLOCK  responds  to  the  real-time  clock  Interrupts  (no. 2) 
at  25  millisecond  Intervals,  reads  anaiog  inputs  through  ADCs  and  per- 
forms the  calculation  required  for  integration  and  analog  output.  At 
Its  end  It  triggers  subroutine  SCHEDULE  (Interrupt  12)  to  determine 
whether  the  run  Is  to  terminate  or  whether  a seeker  sample-hold  update 
is  required.  If  the  latter  Is  required.  Interrupt  8 Is  triggered  to 
activate  subroutine  SEEKER  to  perform  the  sample-hold  operation  on  the 
line  of  sight  rate.  Subroutine  SCHEDULE  also  writes  to  magnetic  tape  a 
record  of  selected  variables  at  time  intervals  determined  by  the  input 
data.  It  should  be  noted  that  all  event  times  within  the  real-time 
section  are  converted  In  routine  STARTRUN  to  an  integral  number  of  25 
millisecond  time  frames.  In  this  way,  the  current  time  Is  maintained 
by  incrementing  a counter  by  I each  time  subroutine  CLOCK  is  entered. 

kl 

Finally,  the  last  part  of  subroutine  SCHEDULE  performs  the  post- 
realtime  operations  of  reading  back  the  magnetic  tape,  printing  the 
results  and  calculating  and  printing  the  intercept  miss  distance. 

A block  diagram  of  the  interlinking  of  the  subroutines  which  com- 
prise the  program  Is  given  In  Figure  4.5.  Subroutine  ABORT  is  triggered 
from  the  analog  console  and  is  used  to  terminate  a sfmu1ation  run 
externally;  ABORT  Is  connected  to  interrupt  number  3.  External 
Interrupt  13  Is  connected  to  an  entry  point  In  the  hybrid  library  named 
FG$RLS  which  has  the  effect,  when  triggered,  of  releasing  foreground 
program  OV  (which  will  usually  be  the  name  under  which  this  program 
Is  executed)  from  the  system. 

A listing  of  each  routine  In  the  program  Is  given  In  Appendix  A 
and  In  Figure  4.6  a flowchart  of  the  program  execution  Is  outlined  to 
Indicate  the  flow  of  control  during  program  execution.  Table  XV 
contains  the  list  of  ADC's  and  OAC's  used  In  the  program  and  their 
variable  definitions. 
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TABLE  XIV  NAMELIST  Input  to  the  Digital  Program 


FORTRAN  I 
NAME 

1 

SYMBOL  ' 

i 

DEFINITION 

DEFAULT 

VALUE 

T 

i TIME  i 

1 

<t 

Initial  value  of  Independent  variable 

0 

1 DT  ! 

6t 

Integration  and  Clock  Interval 

.025  secs 

i TSTOP 

“ 1 

Termination  time  (unless  Intercept  occurs) 

25  secs 

! TACQ 

t 

“ 1 

Time  of  target  acquisition 

o 

, TPSTRT 

l 

t 

Time  to  start  tape  output 

o 

I TAPETIME 

- 

Time  between  tape  output  records 

1 sec 

1 1 PULSE 

- 

Seeker  sample-hold  rate 

20/second 

1 THEMAX 

- 

Seeker  S-H  Input  angle  limit 

2 deg 

PRMAX 

- 

Seeker  S-H  output  angle  limit 

2 deg 

XT 

xT  1 

Target  Initial  X position 

8000  ft 

ZT 

Target  Initial  Z position 

o 

X 

Xm 

Missile  Initial  X position 

0 

I z : 

Zm 

Missile  Initial  Z position 

“500  ft 

VELOCITY 

Vm 

Missile  Initial  total  velocity 

0 

ALFASCAL 

- 

d*  scale  factor 

20° 

DELTSCAL 

- 

6 scale  factor 

20° 

THESCALE 

- 

6r  scale  factor 
Time  of  engine  burnout 

- 

BURNOUT 

tBO 

Mo 

3 secs 

STMASS 

Initial  mass  of  missile 

I*. 6 slugs 

ENDKASS 

m 

Final  mass  (at  burnout) 

3-9  slugs 

STIYY 

lyyo 

Initial  moment  cf  Inertia 

6.1  s .ft~ 

ENDIYY 

lyy 

Final  moment  of  Inertia 

4.9  s.fr 

S 

S 

Aerodynamic  reference  area 

.267  ftZ 

j CBAR 

c 

Aerodynamic  reference  length 

.584  ft 

! QIC 

- 

Initial  value  of  Q, 

0 

j THETA 1C 

- 

Initial  value  of  0£ 

0 

{ THRUST 

^TX  1 

Thrust 

1500  lb. 

GAIN 

1 

? 1 

i 

Gyro  precession  rate  gain 

10/sec 

r 
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Figure  H.5  Subroutine  Linkage  In  the  Program 
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TABLE  XV  OAC  and  ADC  Usage 


DAC  or  .DC 

■ Variable  Definition 

| Scale 
Factor 

t - -* 

; Analog 
Multiplier  i 

i 

DAC  02 

1 qSc2C  n . 

’ 2 1 yy\/  Aerodynamic  pitch  damping 

i 

1 

. 932.8 

i 

o.  ! 

i | 

DAC  03 

i 

qSc 

IyyCm^P,tch  moment  derivative 

373.1 

i 

! ! 

1 

i 

DAC  05 

-a$r 

jyy  Cmg  Vane  moment  derivative 

3.73 

1 

DAC  10 

Wind  axis  rate  term 

0.2 

' ! 

< 

1 

DAC  12 

Wind  axis  rate  term 

0.1 

-1.0 

i 

OAC  22 

”®EIC  an9lfr  initial  value 

180° 

1.0 

1 

DAC  26 

'Olj.  Initial  pitch  rate 

150°/sec 

.1.0 

1 DAC  28 

» 

A*/Ks  Line  of  sight  rate 

2° 

1.0 

! ADC  00 

Cos0£ 

1.0 

1 

ADC  01 

SlnOg 

1.0 

- 

ADC  06 

i 

Angle  of  attack 

.35 

- 

ADC  08 

i 

6 Vane  angle 

1 

.35 

* 

<T>5 
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Figure  ^.6  Program  Overall  Flow  Chart 
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Figure  4.6  Cont. 
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4,4  Analog  Programming  and  Patching  Diagrams 


This  section  contains  the  analog  patching  diagrams  and  potentiometer 
settings  for  the  mathematical  model  sections  assigned  to  the  analog 
computer.  Component  numbers  refer  to  the  Cl -5000  analog  computer. 


TABLE  XVI  Potentiometer  Settings 


Potentiometer 

Number 

Setting 

Potentiometer 

Number 

1 

Setting 

001 

.2000 

' -- 

648 

.5250 

003 

.2151 

049 

.5250 

022 

.8571 

050 

.6852 

025 

.4000 

051 

.7194 

027 

.1237 

068 

.2674 

028 

.2000 

078 

.6000 

029 

.5200 

105 

.1874 

•t 

031 

.8571 

107 

.9550 

036 

.1000 

115 

.6272 

047 

.9999 

129 

.3096 

137 

.5714 

4.$  Results 


The  simulation  described  In  the  preceding  sections  was  checked  out 
and  run  with  the  following  conditions: 

Target  stationary  at  X^.  * 8000  ft,  Zj  - 0. 

Missile  Initial  position  Xm  » 0,  Zm  - -500  ft. 

Hlsslle  Initial  velocity  0,  Initial  0g  • 0,  initial  Q » 0. 

Thrust  ■ 3000  lb,  burnout  at  3 seconds,  and  the  remaining  data  given 
by  the  default  values  In  Table  XIV. 


Results  printed  by  the  digital  program  are  shown  In  Table  XVII  and 
selected  analog  variables  as  functions  of  time,  reproduced  from  strip 
chart  recordings  are  given  In  Figures  4.H  and  4.15* 
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